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ABSTRACT 

The  exact  response  of  thin  shells  loaded  dynamically  is 
extremely  complicated  and  often  cannot  be  expressed  in  closed 
analytical  form.   Approximate  energy  dissipation  solution 
techniques  are  developed  herein  for  rigid-perfectly  plastic 
thin  shells  of  arbitrary  shape  subjected  to  short  duration 
dynamic  loads.   The  effect  of  finite  deflections  is  considered, 
Approximate  results  are  compared  with  some  exact  analytic 
results  for  simple  problems  and  with  some  experimental  results 
for  more  complicated  problems. 
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CHAPTER  1 
INTRODUCTION 

In  recent  years  realization  has  grown  in  the  engineering 
community  that  structural  designs  for  certain  uses  can  be 
improved  by  allowing  some  plastic  deformation. 

Elastic-plastic  analysis  of  structures  can  be  quite 
complex.   It  has  been  found,  however,  that  elastic  deformation 
can  reasonably  be  neglected  if  plastic  deformation  plays  the 
predominant  role  in  a  problem.   Analysis  has  been  simplified 
by  considering  materials  to  be  rigid-plastic. 

Even  with  use  of  the  rigid-plastic  model  of  material 
behavior,  there  are  few  plasticity  problems  for  which  exact 
solutions  have  been  found.   Great  emphasis  has  been  placed 
on  finite  difference  and  finite  element  methods  of  solution. 

This  thesis  investigates  the  dynamic  plastic  behavior 
of  thin,  rigid-perfectly  plastic  shells  and  indicates 
approximate  solution  methods  which  can  be  used  alone  or  as 
an  early  design  aid  in  conjunction  with  more  complete 
numerical  methods. 
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CHAPTER  2 
SOLUTION  TECHNIQUES  FOR  RIGID-PLASTIC  STRUCTURES 

The  stress  level  at  which  plastic  flow  (yielding)  begins 
in  a  material  subjected  to  simple  tension  is  called  the  yield 
stress  or  the  yield  strength,  a..   A  perfectly  plastic 
material  can  suffer  indefinitely  large  strain  once  the  yield 
stress  is  reached.   A  strain  hardening  material  requires 
increased  stress  if  plastic  flow  is  to  continue.   The  yield 
stress  of  a  strain  rate  sensitive  material  depends  upon  the 
rate  at  which  the  material  is  loaded.   The  material  considered 
in  this  thesis  is  rigid-perfectly  plastic. 

If  a  solid  is  subjected  to  a  multiaxial  stress  field, 
the  onset  of  plastic  flov;  depends  on  the  nine  (six 
independent)  stress  tensor  components  a. ..   The  relationship 
of  the  stress  tensor  components  and  the  yield  stress  in  simple 
tension  can  be  expressed  as  a  yield  function 


F(a^.)  =  0  (2.1) 


If  plotted  in  six-dimensional  stress  space,  equation  (2.1) 
is  a  hyper-surface  and  is  called  the  yield  surface. 

Drucker  has  shown  that  the  yield  surface  is  convex  every- 
where and  that  its  shape  governs  the  stress-strain  rate 
relationship  of  the  material  [2-1] . 
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If  the  a. .  axes  are  also  used  as  the  corresponding 
e. .  (plastic  strain  rate)  axes,  the  strain  rate  vector  must 
be  normal  to  the  yield  surface  at  the  point  corresponding  to 
the  state  of  stress.   If  the  stress  point  is  a  corner  or 
vertex  point  of  the  convex  yield  surface,  then  the  associated 
plastic  strain  rate  vector  must  be  a  linear  combination  of 
the  normal  vectors  to  the  surface  at  the  singular  point. 

The  strain  rate  vector  is,  therefore,  proportional  to 
the  gradient  of  the  yield  function  and 


8F{a.  .) 
^ij  =  ^  -^^  ^2.2) 


where  X  is  a  positive  parameter. 

The  principle  of  normality  which  led  to  equation  (2.2) 
is  often  called  the  associated  flow  rule. 

Many  yield  functions  or  yield  criteria  have  been  proposed. 
Two  which  are  commonly  used  and  which  are  in  good  agreement 
with  experimental  results  are  the  VonMises  yield  criterion 
and  the  Tresca  yield  criterion  [2-2], 

The  VonMises  yield  criterion  assumes  that  yielding 
under  multiaxial  stress  conditions  begins  when  the  distortion 
energy  is  the  same  as  that  at  yielding  under  simple  tension. 
The  VonMises  yield  function  in  terms  of  principal  stresses  is 

i[(a^  -  a^)'  +  (a^  -  o ^) ^    +    {o ^   -   a^)^]    -   o^      =   0         (2.3) 
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The  Tresca  yield  criterion  assumes  that  yielding  begins 
when  the  maximum  shear  stress  is  the  same  as  that  at  yield 
under  simple  tension.   The  Tresca  yield  function  can  be 
expressed  as 

Max{|a^  -  G^\,     \o^    -    o^\,     \a^    "C^i  I  )  "  ^n  ^  ^         (2.4) 

If  these  yield  functions  are  plotted  in  Haigh-Westergaard 
stress  space  (principal  stress  coordinate  axes) ,  the  Tresca 
yield  surface  is  found  to  be  a  regular  hexagonal  prism 
inscribed  within  the  VonMises  yield  surface,  a  right  circular 
cylinder  [2-2]  . 

The  state  of  stress  at  a  point  in  a  plate  or  shell  very 
closely  approximates  a  biaxial  state  of  stress  [a ^   =    0) . 
Figure  2.1  shows  the  form  of  the  VonMises  and  Tresca  yield 
curves  in  tv/o-dimensional  stress  space. 

A  perfectly  plastic  material  can  suffer  unlimited 
deformation  once  the  state  of  stress  reaches  any  point  on 
the  yield  surface.   However,  a  structure  made  of  such  a 
material  will  not  suffer  unlimited  deformation  until  enough 
parts  of  the  structure  have  reached  yield  to  produce  a 
collapse  mechanism.  . 

Determination  of  the  exact  static  load  level  which  will 
cause  collapse  of  a  structure  (i.e.,  the  load  at  which  large 
plastic  deformation  just  begins)  is  often  very  difficult  or 
practically  impossible.   The  static  collapse  load  can  be 
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bounded,  however,  by  use  of  the  Upper  Bound  and  the  Lower 
Bound  Theorems  of  Limit  Analysis  [2-3,  2-4,  2-5,  2-6]. 

The  upper  bound  theorem  states  that  a  structure  will 
collapse  under  the  applied  load  if  there  is  any  kinematically 
admissible  pattern  of  plastic  deformation  (i.e.,  one  which 
satisfies  boundary  conditions  and  yields  admissible  strain 
rates)  for  which  the  rate  at  which  the  external  load  does 
work  equals  or  exceeds  the  internal  energy  dissipation  rate. 

The  lower  bound  theorem  states  that  a  structure  will  not 
collapse  (or  will  just  be  at  the  state  of  collapse)  if  an 
equilibrium  set  of  stresses  can  be  found  which  balance  the 
applied  load  and  which  are  everywhere  on  or  within  the  yield 
surface. 

Both  theorems  assume  that  at  the  beginning  of  collapse 
there  is  no  appreciable  change  in  geometry  of  the  structure. 
Results  are  the  same  for  rigid-plastic  material  as  for 
elastic-plastic  material  [2-7] . 

The  methods  of  limit  analysis  have  been  used  successfully 
to  find  or  bound  the  collapse  loads  of  beams  and  framed 
structures  (2-8,  etc.],  circular  plates  [2-9,  2-10,  etc.], 
rectangular  plates  [2-11],  shells  of  revolution  [2-12,  2-13, 
2-14,  etc.],  and  non-symmetric  shells  [2-15]. 

In  the  case  of  plates  and  shells,  generalized  stresses 
and  generalized  strains  (defined  in  Chapter  3)  are  used  in 
place  of  the  usual  stresses  and  strains  defined  at  a  point. 
The  yield  function  and  the  associated  flow  rule  are  then  ex- 
pressed in  terms  of  the  generalized  stresses  and  generalized 

strain  rates. 
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If  the  load  applied  to  a  structure  is  in  excess  of  the 
static  collapse  load,  acceleration  will  result.   Actual 
collapse  of  the  structure  may  be  averted  if  the  duration  of 
loading  is  very  short. 

Exact  solution  of  such  dynamic  problems  is  more  difficult 
than  exact  solution  of  static  collapse  problems.   The  time 
history  of  the  loading  process  may  affect  the  results,  and  the 
shape  of  the  deformation  may  vary  during  the  deformation 
process. 

Few  exact  solutions  have  been  found  for  dynamic  plasticity 
problems  of  plates  and  shells.   Some  solutions  are  exact  only 
for  simplified  yield  criteria. 

Exact  dynamic  solutions  have  been  obtained,  for  example, 
for  circular  plates  [2-16],  cylindrical  shells  [2-17,  2-18], 
and  spherical  shells  [2-19] ,  under  certain  conditions  of 
loading. 

Stable  geometry-change  effects  caused  by  finite 
deflections  have  also  been  incorporated  for  some  problems 
with  beams  [2-20,  2-21],  plates  [2-22,  2-23],  and  shells  of 
revolution  [2-24,  2-25]. 

The  complexity  of  dynamic  problems  in  plasticity  has  led 
to  the  search  for  relatively  simple  approximate  methods  of 
solution.   Techniques  have  been  found  which,  in  principle, 
provide  upper  and  lower  displacement  bounds  for  a  structure 
which  is  loaded  impulsively  [2-26,  2-27].   Unfortunately, 
the  upper  displacement  bound  technique  requires  a  lower 
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bound  (equilibrium)  point  load  collapse  solution,  which  is 
extremely  difficult  to  find  for  complex  structures. 

An  approximate  method  for  describing  the  dynamic 
behavior  of  rigid-plastic  beams  and  flat  plates  was  recently 
developed  by  Jones  [2-28] .   The  method  assumes  a  time- 
independent  deformation  mode  shape,  ignoring  the  motion  of 
travelling  hinges,  and  accounts  for  the  effects  of  finite 
transverse  deflections.   Comparison  with  experimental  results 
in  [2-28]  and  in  [2-23]  showed  good  corelation. 

The  major  development  of  this  thesis  is  an  extension  of 
Jones'  approximate  method  for  beams  and  plates  to  the 
dynamic  behavior  of  thin,  rigid-plastic  shells. 
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CHAPTER  3 
BASIC  EQUATIONS  FOR  AN  ARBITPJ^RY  THIN  SHELL 

Section  1.   Assumptions 

The  development  of  certain  sets  of  equations  basic  to 
the  deformation  process  of  an  arbitrary  shaped  shell  is  out- 
lined in  this  chapter.   The  equations  developed  are  those 
which  relate  generalized  strain  rates  to  the  displacement 
rates  of  the  shell  (strain/curvature  rate  equations)  and 
those  which  relate  the  generalized  stresses  and  their  deri- 
vatives to  the  external  load  and  the  D'Alembert  forces  (the 
equilibrium  equations  or  equations  of  motion) . 

Several  simplifying  assumptions  have  been  made  in  the 
development.   The  results  are  applicable  to  the  wide  class 
of  problems  for  which  the  transverse  deformation  (perpendi- 
cular to  the  mid  surface  of  a  shell)  is  the  predominant 
deformation. 

The  major  assumptions  are  listed  below. 

1.  The  shell  is  thin;  i.e.,  the  smallest  radius  to 
wall  thickness  ratio  is  greater  than  or  equal  to 
about  ten. 

2.  Material  is  homogeneous  and  isotropic. 

3.  Straight  lines  perpendicular  to  the  undeformed 
midsurface  of  the  shell  remain  straight  and 
perpendicular  to  the  deformed  midsurface.   This 
is  usually  known  as  the  Love-Kirchof f  assumption. 
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4.  Displacements  in  the  transverse  direction  and 
their  derivatives  are  much  larger  than  the 
displacements  and  derivatives  of  displacements 
tangential  to  the  shell  midsurface. 

5.  Transverse  displacements  may  be  moderately  large, 
but  the  strains  which  they  produce  are  infinitesimal. 

6.  The  effects  of  rotary  inertia  are  small  compared  to 
those  of  linear  or  translational  inertia. 

7.  The  normal  stress  in  the  transverse  direction  is 
small  and  is  assumed  to  be  zero. 

Section  2.   Differential  Geometry 

In  order  to  find  the  strain  rate,  curvature  rate,  and 
equilibrium  equations  one  must  be  able  to  describe  the 
geometry  of  the  shell.   Therefore,  some  space  will  be  devoted 
to  the  differential  geometry  of  a  curved  surface  (the  mid- 
surface  of  a  shell,  for  example). 

The  midsurface  of  a  shell  is  a  two-dimensional  space. 
Some  information  about  the  surface  is  intrinsic,  that  is,  it 
can  be  expressed  in  terms  of  the  two-dimensional  surface 
alone.   For  a  complete  representation  of  the  surface,  however, 
other  information  is  required  which  relates  the  surface  to 
the  three-dimensional  space.  [3-1] 

Consider  the  surface  element  shown  in  Figure  (3.1). 
Each  point  in  the  surface  can  be  represented  by  a  three- 
dimensional  position  vector  r  (with  components  r  )  with 
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respect  to  the  rectangular  cartesian  coordinates  X  .   Each 
point  can  also  be  represented  by  its  coordinates  with  respect 
to  a  curvilinear  coordinate  system  which  lies  in  the  surface. 
The  most  simple  surface  coordinate  system,  x]    ,    is  the  family 
of  mutually  orthogonal  surface  curves  which  correspond  to  the 

curves  of  maximum  and  minimum  normal  curvature.   The  spatial 

i     i   ct 
coordinates  can  be  expressed  as  X   =  X  (ri  )  , 

It  should  be  noted  that  Roman  indices  (i,  j,  k,  etc.) 

will  be  used  to  denote  spatial  tensor  variance  and  will  have 

the.  range  1  to  3 .   Surface  variance  will  be  denoted  by  Greek": 

indices  {a,^,y,    etc.)  which  will  have  the  range  1  to  2. 

Consider  two  adjacent  points,  P,  and  P„,  lying  in  the 

ct  . 

surface  defined  by  ri  .   As  the  two  points  approach  each  other 

closely,  the  length  of  the  spatial  differential  vector  dr 

approaches  the  differential  curve  element  length  ds . 

Thus  in  the  limit  as  ds  approaches  zero, 

(ds)^  =  dr  •  dr  (3.1) 

By  use  of  the  Chain  Rule  of  differentiation. 


and 


dr   =  r,   dri  (3.2) 

a 


dr  •  dr  =  r^,   r^,^  dn^'dn^        •  (3.3) 

a    3 


Therefore 


(ds)2  =  r^,   r^,,  dn°^dn^  (3.4a) 

a    3 
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or 


(ds)2  =  a  ^dn°^dn^  (3.4b) 


where 


a  ^  =  r  ,   r  / ^   =  surface  metric  tensor 
a3,      a    3 


Since  the  n   are  orthogonal,  a^ ^  =  a^,  =  0.   Equation 
(3.4b)  is  called  the  first  fundamental  form  of  the  surface. 

The  surface  coordinate  system  also  has  a  set  of  base 
vectors,  a  .   These  base  vectors  are  tangent  to  the  surface 
m  the  direction  of  the  corresponding  surface  coordinate  x]    . 

The  surface  base  vectors  are  defined  as 


a      =   ^^  (3.5) 


The  surface  base  vectors  are  needed  to  relate  the  surface 
to  three-dimensional  space.   This  is  done  by  forming  the  unit 
normal  vector  to  the  surface.   Since  the  a   are  tangent  to 
the  surface,  their  cross  product  will  be  normal  to  the  surface 
and 

a   X  a_ 

a     3 


n  =  I-    -  ,  (3.6) 

a   X  a„ 

'  a    3 


with  the  order  of  multiplication  chosen  so  that,  by  convention, 
the  direction  of  n  is  from  the  concave  side  toward  the  convex 
side. 
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The  additional  relationship  needed  to  represent  a  surface 
is  that  which  relates  the  change  in  the  nonnal  vector,  n,  to 
change  of  the  position  vector,  r,  as  one  progresses  along 
the  surface.   This  relationship  is  called  the  second  funda- 
mental form  of  the  surface  and  is  expressed  as 


dn  .  dr  =  b  ^  dn^'dn^  (3.7) 

ab 


where 


-  1  .Jn       "5r    "5n    "?r  , 


8n    5r)     3n     8n 

In  order  to  find  the  elements  b  „,  one  must  differentiate 
the  unit  normal  vector  n  with  respect  to  the  curvilinear  surface 
coordinates  ri  . 

Recall  that  a   = 


9n 


and  that  r  is  the  position  vector 


r  =  b.r^  =  b.X^ 
1       1 


Since  the  spatial  base  vectors,  b.,  for  a  rectangular  cartesian 
system  are  constants. 


where 


X-'  denotes 


3X^ 


a         ^  a 
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which  is  the  rate  of  change  of  the  spatial  coordinate  X^  with 
respect  to  that  of  the  surface  coordinate  t]^ , 

The  numerator  of  equation  (3.6)  using  the  orientation  of 
Figure  (3.1)  is 


a^^  X  32  =  (b.  X^)  X  (bj^  X^)  (3.9a) 


which  is  expressible  in  tensor  notation  as 


^1  '^  ^2  =  ^ijk  4   4  '^i  '3.9b) 


where  £...,     is  the  alternator  symbol. 
Ilk  ^ 


e. .,  =1       for  ijk  an  even  permutation  of  123 

£. .,  =  -1       for  ijk  an  odd  permutation  of  123 

e.  .,  =0       for  two  or  more  indices  being  equal. 

The  denominator  of  equation  (3.6)  is  the  magnitude  of  the 
numerator. 


a.,  x  a„  I  =  |a,|  [a^j  sin9  (3.10) 


OL 

Since  the  n   are  specified  to  be  mutually  orthogonal, 
sin9  =  1. 

Furthermore,  since  the  magnitude  of  a  vector  is  the  square 
root  of  the  scalar  product  of  the  vector  with  itself. 
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|a  I  =  [r^-,   r"-,  ]^  =  /T" 
'a'        a     a       aa 

(No  summation  on  a) 

Equation  (3.10)  therefore  reduces  to 

I-     -  I     r         A 

'^1  ^  ^2'  "^  ^^11  ^22^  (3.11) 

and  equation  (3.6)  takes  the  form 

e.  .  •  X^  X^  b. 

7:   -      Ilk   1   2   1  /-)  1  n\ 

n  =  ' = 1 (3.12) 

^^11  ^22^^ 

To  find  the  absolute  rate  of  change  of  the  normal  vector 

with  respect  to  the  change  of  position  on  the  curved  surface 

one  must  subtract  the  apparent  change  due  to  the  curvature  of 

ct 
the  n   coordinate  system.   Thus  the  tensor  (or  covariant) 

derivatives  of  n  must  be  used. 


n/a  =  b.  n  /a  (3.13) 


where 

n  /a  =  n  ,   +  r  .,  n-*  X 
a     jk     a 

and  r .,  is  the  Christoffel  symbol  of  the  second  kind  for  the 
three-dimensional  space  X  . 

The  X   being  used  form  a  rectangular  cartesian  coordinate 
system  and 

r\  =  0. 


26 


Therefore,  the  covariant  derivative  of  a  space  vector 
with  respect  to  the  surface  coordinates  is  identically  equal 
to  the  partial  derivative  of  the  space  vector  with  respect  to 
the  surface  coordinates. 

That  is 


n  /a  =  n  ,  (3.14) 


Equation  (3.7)  can  be  expanded  using  (3.8)  as 


dn  •  dr  =  (n/a  dn^')  •  (b .  X^  dn^)  (3.15a) 

3   P 


which  in  turn,  by  using  equation  (3.14),  may  be  cast  in  the 
form 


dn  •  dr  =  n^,   X^„  dn"dn^  (3.15b) 

a    p 


Comparing  equations  (3.15b)  and  (3.7)  shows  that  b  g  may 
be  found  from 


b„g  =  n\^   xj  (3.16) 


The  normal  curvature  of  a  surface  curve  is  related  to 

the  tensors  a  ,  and  b  ^  [3-1] .   Passing  through  every  point 
ab      a  3 

on  the  surface  there  are  two  mutually  orthogonal  curve  segments 
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whose  normal  curvatures  have  extreme  values.   If  the  coordinate 

ct 
curves  n   are  aligned  with  these  lines  of  principal  curvature. 


''la  =  ''21  =  0 


and  the  principal  curvatures  are 


1     ^11  1     ^22 

|-  =  -^     and     I"  =  z^  (3.17) 

^1    ^11  ^2    ^22 


Where  R  is  the  radius  of  curvature  of  the  coordinate 
a 

a 
curve  T]    . 

In  order  to  fully  describe  a  shell  one  must  bring  in  a 
third  local  coordinate  normal  to  the  midsurface.   This 
direction,  Z  or  n ^ ^  is  used  to  express  shell  thickness. 

Therefore  we  shall  now  consider  a  family  of  surfaces 

which  are  parallel  to  the  shell  midsurface.   That  is,  surfaces 
each  of  which  is  a  constant  distance  Z  away  from  the  midsurface 

The  sense  of  Z  agrees  with  that  of  the  unit  normal  vector  n. 

The  surface  metric  tensor  of  a  parallel  surface  is 

designated  a  _.   Since  the  surface  is  parallel  to  the  raid- 
ex  p 

surface,  it  has  the  same  principal  directions  and 


^  /\  ** 

^12  =  ^21  =  °- 


The  nonzero  values  of  a  „  are  found  from  the  application 

dtp 

of  equation  (3.1)  to  the  parallel  surface  at  distance  Z  from 
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the  midsurface.   The  position  vector  f  to  a  point  on  the 
parallel  surface  is 


r^  =  r^  +  Zn^  (3.18) 


and 


df^  =  (r^,^  +  Zn^,^)  dn^'  (3.19) 


Proceeding  as  in  equations  (3.3)  and  (3.4)  one  defines 


h&   =  '-^^'a  -^  ^""'a'  <^^'6  ^  '"  'b'  <^-2°^' 


which,  when  the  indicated  multiplication  has  been  carried  out, 
becomes 


a  o  =  r  ,   r  ,  o  +  Zr  ,   n  ^^  +  Zr  ,„  n  ,   +  Z  n  ,   n  ,o 
a3      a    3       a    3       3    a        a    3 

(3.20b) 

Since  the  position  vector  component  r   is  the  coordinate  X  , 
the  second  and  third  terms  of  equation  (3.20b)  may  be  expressed 
through  the  use  of  equation  (3.16)  as 


Zr^,   n^,^  =  Zb^  (3.21) 
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The  Weingarten  formula  [3-1]  for  a  unit  normal  vector 
directed  outward  states  that 


n^,   =  a*^^  b;,   X^  (3.22a) 

a        oa   Y 


which  permits  expression  of  the  third  term  of  equation  (3.20b) 
as 

i     i       Sy  ,    ^i  Xy  ,    ^i     <5y   Xy  ,    , 

n,   n,_  =  a'b.   X  a^  b,  „  X   =  a  '  a   b,.   b,  „  a 

a    3        6a   Y  X3   y            (5a   X3  y\i 

(3.22b) 
since 


X^  X^  =  a 
Y  y    Yy 


a 
For  the  principal,  orthogonal  coordinate  system  r\    , 


^up  ^  _x_  for  a  =  B 

^aS 
and 

3°^^  =  a  ^   =  0  for  a  7^  B. 

a3 

Therefore 

b^ 

n^,   n^,   =  — ^  (no  sum  on  a)  (3.23) 

a     a    a 

aa 

Applying  equations  (3.21)  and  (3.23)  to  equation  (3.20b) 
yields  the  result 
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b  b 

aa        aa 
(no  sum  on  a)  (3.24a) 


However 


b 

aa 


a      R 

aa    a 


(no  sum) 


and,  therefore. 


^aa  =  ^aa^^  ""  h^  ^  ^^°  ^^  °^  °^^  (3.24b) 


With  this  result  a  shell  of  finite  thickness  can  be 
described  in  terms  of  the  first  and  second  quadratic  forms 
(and  tensors)  of  the  shell  middle  surface  and  directed 
distance  normal  to  the  midsurface. 

Section  3.   Strain  Rate  and  Curvature  Rate 

The  geometric  representation  of  a  shell  developed  in  the 
preceding  section  may  now  be  combined  with  the  definition  of 
strain.   The  approximations  of  Section  1  are  then  used  to 
develop  the  generalized  strains  of  a  thin  shell. 

Consider  the  local  normal  direction,  the  Z  coordinate. 
It  is  also  the  third  coordinate  in  the  right-handed,  orthogonal 
curvilinear  system  n\  r]^  ,    n  ^   =  Z.   The  corresponding  metric 
tensor  component  is  a^.,  =  1. 
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As  the  shell  deforms,  let  the  coordinate  system  also 
deform  such  that  the  coordinates  after  deformation  are 
numerically  equal  to  the  coordinates  before  deformation.   In 
general,  the  metric  tensor  will  change  to  accommodate  the 
deformation.   Since  the  normals  to  the  midsurface  remain 
normal  during  deformation,  ^-.o'  ^-31'  ^0-3  ^^*^  ^-50  ^iH 
remain  zero. 

Interest  has  been  limited  to  finite-deflection  with 
infinitesimal  strains,  and  for  convenience  the  undeformed 
shape  is  used  as  reference.   The  description  of  motion  is 
therefore  Lagrangian. 

The  Lagrangian  strain  tensor  in  terms  of  the  three- 
dimensional  curvilinear  coordinate  system  n   and  in  terms 
of  the  covariant  and  contravariant  vector  displacement 
components  y .  and  y   is  [3-2]: 

E.  .  =  i(y ./.  +  VI.  /^  +  y^^^  y,  /J  (3.25) 


:..=7(y./-+y-/-+y    y,/- 


where  (   )/„  denotes  covariant  differentiation  of  (   )  with 
respect  to  n  . 

Since  the  local  three-dimensional  coordinate  system 
n   is  orthogonal, 

a^^  =  J-  (3.26a) 


and 


a. 
11 


a.  .y^  =  y.  (3.26b) 

11^      1 


(no  sum  on  i) 
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The  form  of  the  covariant  differentiations  referred  to 
in  equation  3.25  is 


^j/i  =  ^j,i  -  r-i  U,  (3.27a) 


and 

y^/i  =  y^,  +  tI    y^ 


£i  y  (3.27b) 


where  (  ),.    indicates    . '    and  the  Christoffel  symbol  T        is 
1  ti^i  mn 

on 


ii  £ 

But  a  -"  =  a.  .  =0  for  i  7^  j  and  T        is  different  from  zero 
13  -^      mn 

only  if  k  =  £.   Then 


r^   =  ■^a^^{a  ,    +  a  „    -  a    „}  (3.28b) 

mn    ^      m£,n    n£,m    mn,£  ' 

(no  sum  on  £) 


Furthermore,  if  Z,    m,    and  n  are  all  different,  all  three 

terms  in  the  bracket  of  equation  (3.28b)  are  equal  to  zero 

and  thus  r    is  zero. 
mn 

Z 
Therefore  T         is   different  from  zero  if  and  only  if  at 

mn  -^ 

least  two  of  the  indices  Z,    m,  and  n  are  equal. 

The  coordinate  curve  n^  =  Z  is  a  straight  line  and 

a--  =  1.   So  a-,-,  .  =  0.   Also  since  a . -,  =  0  for  i  5^  3,. 
33  33 ,;]  i3  ■  ' 

a.  -.  .  =  0. 
i3,:i 
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0 

The  three  components  of  T        can  therefore  be  written  as 

mn 


mn  ml,n  nl/in  mn,l 


•p2         1^22/^  j.^  '^  \ 

^mn   "   ^    ^       ^^m2,n   ^    ^n2,m   "    ^mn,2^ 


r^      =   -i    a        .  (3.29) 

mn  ^      mn,3  ^  ' 


Let  us  look  now  at  one  term  of  the  strain  tensor  E. .. 

ID 

^11  =  ^1/1  ^  ^^^"^"^^   ^1/1  ^  ^'^^  ^2/1  ^  ^^''^  ^3/1^   ^^-^0) 


Before  expanding  this  representative  term,  let  us  invoke 
the  assumption  that  transverse  deflections  w  are  much  larger 
than  the  in-surface  physical  deflections  u^  and  u^.   There- 
fore, the  only  terms  retained  v/hich  involve  products  of 
displacements  (or  products  which  indicate  derivatives  of 
displacements)  are  those  for  which  w  is  involved  in  each 
factor.   Because  the  Z  coordinate  is  a  straight  line, 
w^  =  w-  =  w.   With  this  fact  and  the  approximation  noted 
above,  rewrite  equation  (3.30) as 


^11  "  ^1/1  "^  ^  ^/l  ^/l  (3.31a) 


which  becomes 


^k        .   1  /_    X  2 


^11  "  ^1  1  ~  ^11  ^k  ^  7(w^^)^  (3.31b) 
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after  performing  the  covariant  differentiation  and 
eliminating  the  mixed  product  displacement  (because  of  the 
assumption  that  transverse  displacements  are  much  larger  than 
in-surface  displacements) . 

Expanding  this  by  using  equation  (3.29)  yields 

^11  -  ^1,1    2 ""  2 ^ 

hi    3      (w  t)^ 
^i'w   +      — ^  (3.31c) 


The  relationship  between  the  strain  tensor  E. .  and  the 
physical  strain  tensor  e. .  is  derived  from  the  expression  [3-2] 


2E.  . 

^ 1  =  (1+e.  . )  (1+e  .  .)  sin  y.  .  -  1  (3.3  2a) 

[a.,  a..]^       ^^      3D       'id 

11   DD 

(no  sum) 
where  6. .  is  the  amount  by  which  the  right  angle  between 
coordinate  curves  n   and  r\-^    changed  during  deformation. 

For  infinitesimal  strains,  e,,  <<  1     (no  sum)   and 
the  physical  strain  tensor  can  be  expressed  as 


E.  . 

£.  .    =  ^ ^         (no  sum)  (3.32b) 

^3    [a.  .  a.  .]* 
11  33 

where  e  .  .  +  e  .  .    =   y  .  .  for  i  5^  j  . 
ID     Di     ID 
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The  vector  displacements  y  are  related  to  the  physical 
displacements  by  [3-1] 


u.    =/7^  y.    =   -^ —  (3.33a) 

1  1  rr 

/a.  . 
11 

and 


u.    =    v/177      y^    =        y^  (3.33b) 

1  11      "^  -^ 


11 


(no   sum  on   i) 
andj.    as    noted   earlier. 


w   =   w^    =  w-  (3. 33c) 


The  physical  strain  tensor  component  e, .  is  therefore 


_  ^-^  ^l\l    ^11.1  ^1   ,  ^11,2  ^2   , 
^11 ^  ^ 


^11        ^  ^11*^^11    ^^^22  ^11 


^11,3  ^   ^  i\li' 

^  ^11  ^  ^11 


The  remaining  physical  strain  tensor  components  are  found  in 
the  same  manner.  The  nine  (six  independent)  physical  strain 
tensor  components  are 


36 


11 


2*^  ^11 


2^^  ^11 


^11.3  ^  ^  ^^.1^' 
2  a 


11 


2  a 


11 


(3.34a) 


22 


^^"^  ^2^  ,2 


'22 


^22,2  ^2   ^  ^22,1  ^1   _j_ 


2/a^  ^22 


"^"^^11  ^22 


^21^  w+  (w^2^^ 


2  a 


22 


2/\ 
a 


22 


(3.34b) 


C33  =  0 


(3.34c) 


^12    ^21 


2*^^11  ^22 


^  5        n 

11   P  ??  1    2 


/a 


11 


/I 


22 


^^^11  ^22 


^•^^11  ^22 


(3.34d) 
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(a   u  )      w       a     u 

e,-,  =  £-,,  =  — — — i— ^  +  -^ ^^^^ — -  (3.34e) 

2/a,,        ^*^^ii        11 


^"^^  "^2^  3    ^2      ^22  3  ^2 
£-_  =  £--  =  ^ ^—^  +   -^^ ^^-^-^^ ^         (3.34f) 

^  J  J  ^  —  /A  —  /A  ^ 

2/322        2*^^22       ^22 


Further  simplifications  can  be  made  through  use  of  the 
Love-Kirchof f  assumption  that  straight  lines  normal  to  the 
undeformed  midsurface  remain  straight  and  normal  to  the 
deformed  midsurface.   The  in-surface  components  of  displace- 
ment u   can  be  expressed  in  terms  of  the  displacement  u  r.  at 
the  midsurface,  the  slope  (u   Z) ^  of  the  rotated  normal,  and 
the  directed  distance  Z  of  the  point  along  the  normal.   That 
is 

Recall  that  the  surface  metric  tensor  at  height  Z  is 
related  to  the  midsurface  metric  tensor  by  equation  (3.24b). 


Si        =  a   (1  +  |-)2  (3.24b) 

a 


(no  sum  on  a) 

The  rate  of  change  of  this  tensor  element  with  respect 
to  Z  is 

\a,Z  =  2a^^(l  +  |-)  |-       (no  sum  on  a)  (3.36) 

a    a 
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The  limitation  of  thin  shells  stipulates  that 


g 


and,  therefore /  that 


g 


Equations  (3.24b)  and  (3.36)  are  simplified  to 


a        =  a 

gg          gg 

and 

2a 

a            _        gg 

ga,Z           R 
a 

(3.37a) 


(3.37b) 


Since  the  normals  to  the  midsurface  remain  normal, 


^13    ^23    °* 


By  carrying  out  the  indicated  operations  of  equation  (3.34e) 
with  the  approximations  of  equation  (3.37)  and  by  evaluating 
the  results  at  the  midsurface  (Z=0) /  one  finds  that 


^  1    "lO 
(-1,Z^0  =  -  7^+^  (3.38a) 


and 


"^  2     "^20 
(-2,Z^0  =  -  7:l=^-#  (3.38b) 

/a22      2 
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Substitution  of  equations  (3.38)  into  equation  (3.35) 

yields 

u     w 

u,  =  u  (,  +  Z(-i2. d^)  (3.39a) 

1    ,^ 

and 

""20    "^  2 
-2  =  -20  ^  ^(nf  -  7=^  (^-^^^^ 

2     /^ 

Let  us  now  substitute  equations  (3.39)  along  with  equations 
(3.37)  into  the  expression  for  e-,w  equation  (3.34a), 
collecting  those  terms  of  which  Z  is  a  factor. 


=  ^^^^l0^1  _  ^11.1  ^10  ^  ^11.2  -20  ^ 
^11        2/^  a^^    2  a^3_/I^ 


/a   u 
(w   )  2       ( — ^ ^) 

I-  + — ^^1-]  +  ,,__^i__:i_    ^,11    ^11.1  ^10  , 

^         "^^11  ^11  ^11     2(a^^)^  R^ 


^11.1  ^,1  ,  ^11.2  -20      ^11.2  ^.2  ^3^4Qj 

2(a^^)2      2a3_3_/i^  R2    2a^^  a^^ 


This  result  and  the  similar  results  obtained  from 

equations  (3.34  b  and  d)  will  be  expressed  in  terms  of  the 

commonly  used  generalized  strains  e  „  the  midsurface  normal 
-^  ^  otO 

strains,  Yn  the  midsurface  shear  strain,  k      the  bending 
0  a  ^ 

changes  of  shell  curvature,  and  k-„  the  twisting  change  in 
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shell  curvature,  which  are  defined  by  the  following  set  of 
equations. 


^11  "^  ^10  "*"  ^"^1  (3.41a) 


^22  =  ^20  ■"  ^^2  (3.41b) 


^12  ■"  ^21  =  "^0  -^  2ZK^2  (3.41c) 

Before  writing  the  three  strain  expressions  and  the  three 

curvature  expressions,  let  us  define  two  metric  scale  factors 
A  and  B  such  that 

A(n\n')  =  [a^^(n\n')]^  (3.42a) 


and 


B(n\n')  =  [a^^{^' ,r]^)]^  (3.42b) 


Proceeding  as  in  equation  (3.40)  and  using  the  notation  of 

equations  (3.42),  the  (midsurface)  strain  components  and 

(change  of)  curvature  components  defined  by  equations  (3.41) 

are  found  to  be 

u,       A   u^  (w  ^  )  2 

^10 A—^  ^  ~AB        R^     2A^  (3.43a) 

u^-  ^    B  -  u- -  (w  ^) ^ 
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,  =l,l^.^l,    .^l,!lo.^l,  ,3.43e) 

-)       A  ,  10     ,1.    .  B  ,  20     ,2.  ,->  /i-)^\ 

2^12  =  B  ^AR^  -  -A^^2  -^  A  ^BR^  -  "B^^  ,1  ^^'^^^^ 


It  should  be  noted  that  retention  of  finite  (moderately 
large)  transverse  deflections  gives  rise  to  the  nonlinear 
term  in e^ - ,  e^.  and  y ^.      These  nonlinear  terms  are  missing 
if  all  deflections  are  infinitesimal. 

The  strain  rate  and  curvature  rate  equations  are 
obtained  by  differentiating  each  of  equations  (3.43)  with 
respect  to  time.   In  the  case  of  infinitesimal  deflections, 
the  strain  and  curvature  equations  are  homogeneous  with 
respect  to  time.   If  finite  deflections  are  accounted  for, 
the  nonlinear  term  in  each  of  the  midsurface  strains  produces 
a  nonlinear  term  in  the  corresponding  strain  rate  equation. 

The  strain  rate  and  curvature  rate  equations  are 


u,  _  T  A  ^  u_-  •     w  T  w  T 

^10      A        AB  ^1    ^  IJ.^^a; 

•  •  • 

u  B    u  •WW 
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^0         B^    A    \2         A^    B    \l  AB  AB  IJ.^'ic; 


-,    u  w  A        u  w 


nU^rt  W^  B-U--  w_ 

<2   =  l'-|f-^',2-KI^<Hf     -^'  <3-44e) 


2^12    =    B^AR^   -    A^\2    "■   A^BRJ   "    B^^  1  ^^'^^^^ 


The  fact  that  the  strain  and  curvature  expressions 
(and,  therefore,  the  rate  equations)  derived  in  this  chapter 
contain  approximations  should  be  reemphasized. 

There  have  been  many  attempts  to  derive  strain  and  cur- 
vature expressions  for  thin  shells  using  the  Love-Kirchof f 
assumption  or  something  equivalent,  and  there  are  almost  as 
many  different  results  as  there  have  been  attempts.   Differences 
creep  into  the  curvature  expressions  due  to  the  approximations 
made  and  due  to  the  stage  at  which  the  approximations  are 
made.   The  author  has  not  found  any  disagreement  in  the 
expressions  for  (midsurface)  strain. 

Table  3.1  gives  a  comparison  between  some  representative 

theories.   The  Table  shows  what  terms,  if  any,  must  be  added 

to  the  corresponding  curvature  expression  of  this  thesis  to 

bring  it  into  agreement  with  the  theory  referenced.   Symbols 

and  directions  have  been  changed  to  agree  with  those  of  this 

thesis. 
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Budiansky  and  Sanders  [3-8]  defined  a  pseudo-curvature 
in  order  to  produce  a  symmetric  tensor  from  which  curvature 
could  be  found.   The  comparison  in  Table  3.1  was  made  with 
the  curvature,  not  the  pseudo-curvature. 

The  various  curvature  results  would  lead  to  differing 
results  for  the  equilibrium  equations  as  well.   Only  the 
set  of  equilibrium  equations  which  is  consistent  with  strain 
rate  and  curvature  rate  equations  of  this  thesis  will  be 
presented.   The  finite  deflections  considered  do  not  affect 
the  curvature  rate  expressions. 
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Section  4.   Equilibrium  Equations 

Let  us  now  investigate  the  equilibrium  of  a  differential 
element  of  a  thin  shell.   In  order  to  ensure  that  the 
equilibrium  equation  results  will  have  approximations 
consistent  with  the  rate  equations  (3.44)/  the  equilibrium 
equations  are  developed  from  an  energy  approach  rather  than 
from  a  strictly  geometric  approach  [3-10] . 

In  the  theory  of  plates  and  shells  the  resultant  forces 
and  moments  which  act  on  the  midsurface  at  an  element  face 
are  normally  used  instead  of  the  stresses  themselves.   The 
resultants  are  defined  by  a  force  or  moment  balance  and  are 
called  generalized  stresses. 

A  typical  shell  differential  element  is  shown  in 
Figure  (3.2)  and  indicates  the  positive  sense  of  stresses, 
forces,  moments,  and  pressures  used  in  this  thesis. 

Taking  a  force  balance  in  the  n ^  direction  on  the  face 
perpendicular  to  n ^  yields 


N-Bdn'  =   /^^^  a,  Bdn'  -i dZ 

^  -h/2    ^  ^2 


A  moment  balance  in  the  t\^    direction  on  the  n  ^  face  yields 

R  +7 

M^Bdn^  =  /^^^  o.   Bdn^  z  -^ —  dz 

^  -h/2    ^  ^2 
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Fig.    3.2 
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And  a  force  balance  in  the  Z  direction  on  the  n ^  face  yields 


R  +Z 

Q.Bdn'  =  /^/^  T,^  Bdn'  -|—  dz 

^         -h/2   ^-^         ^2 


Applying  the  same  technique  to  the  remaining  forces 
and  moments  and  using  the  thin  shell  assumption  that  —  <<  1, 
the  definitions  of  the  generalized  stresses  become 


N   =  /^^■^   a^dZ  (3.45a) 

^        -h/2  ^ 


N^  =  /  a2dZ  (3.45b) 


^12  "  ^21  "  ^      "^12^^  (3.45c) 


M  =  /    a  ZdZ  (3.45d) 


M  =  /   a2ZdZ  (3.45e) 


^12  "^  ~^21  "^  •''    "^12^^^  (3.45f) 


Q^  =   /  -^12^^  (3.45g) 


Q2  =   /   "^23^^  (3.45h) 


The  integration  limits  are  the  same  in  each  definition, 
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Because  of  the  Love-Kirchof f  assumption  that  normals  to 
the  midsurface  remain  normal  and  straight,  e,^  =  e^^  =  0. 
Therefore  there  is  no  transverse  shear  strain  energy,  and  Q, 
and  Q-  are  considered  to  be  reaction  forces  necessary  for 
transverse  equilibrium. 

The  affect  of  the  shell  thickness  has  been  taken  into 
account  in  the  definitions  of  equation  (3.45)  and  the 
resultant  forces  and  moments  are  assumed  to  act  on  the  shell 
midsurface. 

Figure  (3.3)  illustrates  the  variation  of  the  generalized 
stresses  and  reactions  of  a  loaded  shell  element. 

The  principal  of  virtual  velocities  states  that  if  any 
stress  field  which  is  in  equilibrium  with  the  applied  loads 
acts  in  conjunction  with  any  arbitrary,  continuous  virtual 
velocity  field  which  satisfies  the  actual  boundary  conditions, 
then  the  external  work  rate  (due  to  reaction  forces  and  applied 
load)  is  equal  to  the  internal  energy  dissipation  rate. 

Symbolically  this  is  represented  as 


W   ^      =  W.  ^     T  (3.46) 

external    internal 


The  internal  energy  dissipation  rate  is  determined  by 
integrating  the  strain  energy  rate  over  the  volume  of  the 
body.   This  result  can  also  be  obtained  by  integrating  the 
product  of  the  generalized  stresses  with  the  corresponding 
generalized  strain  rate  over  the  undeformed  midsurface  area, 
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Fig.    3.3 
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^int  =  ^f^l^l  ^  ^2^2  -^  ^12^0  ^  ^1^1  -^  ^2^2  ^    2M^2^12^^S 

(3.47a) 


The  external  work  consists  of  contributions  from  the 
work  done  by  boundary  forces  and  moments,  by  the  D'Alembert 
forces  (inertial  forces) ,  and  by  the  load. 

Figure  (3.4)  indicates  the  boundary  forces,  moments 
and  associated  displacement.   The  boundary  is  a  curve  and 
the  boundary  work  rate  is,  therefore  a  line  integral. 


ext      aB   a3    a  a    a6  a3 
boundary 


/P.u.  dA  -  /yii.u.dS  (3.47b) 

S  ^  ^       S   ^  ^ 


First  let  us  evaluate  a  representative  term  of 
equation  (3.47a)  by  substituting  equations  (3.44)  for  the 
strain  and  curvature  rates. 

The  first  term  of  equation  (3.47a)  becomes 


s         n  n  1 


(Adn^Bdn^)  (3.48) 
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Fig.    3.4 
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The  first  and  last  terms  on  the  right  hand  side  of 

equation  (3.47)  contain  surface  derivatives  of  displacement 

rates.   Let  us  integrate  those  two  terms  by  parts  in  order 

to  produce  a  line  integral  at  the  boundary  which  can  match 

that  found  in  W   , .* 

ext 

N^u^  «  ,    N^w  ,w 


//[  V°'^  "^   \^^  ^^]  ABdn^dn2  = 


N - Bw  - w  ^ 


N   Bw 
/{[N^Bu^q]    +    [-A__lAw]}    dn2    - 


N   Bw    ^     • 

where  the  expression  within  the  {}  is  to  be  evaluated  at  the 
shell  boundary.   The  integral  involving  that  expression  can 
be  expressed  as  the  line  integral 


N , Bw  , w 
boundary 


The  first  term  of  equation  (3.47)  can  therefore  be  written 
as 


*For  ease  in  presentation  the  n   symbols  will  be  typed  as  n 
and  the  integration  paths  will  be  omitted  from  the  //. 
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N  -  Bw  ^  W 


N^ABw 


N  Bw  ^ 
(  \  '    )  ^^   w]  drij^  dn2  (3.49) 


Next  let  us  consider  the  components  of  motion  at  a 
shell  boundary.   Figure  (3.5)  shows  a  differential  profile 
of  a  shell  midsurface  at  a  boundary  before  and  after  defor- 
mation and  indicates  the  displacement  rate  components  on  the 
edge  ri-,  =  constant. 

The  angle  <])-  is  the  sum  of  the  angles  a  and  b. 

^10    "^  1 
*ll  =  -i^-^  ^'-^'-^ 

The  corresponding  angle  for  the  edge  n-,  =  constant  is 

"^20    "^  2 
*22  =-|f  -  -^  (^-^^^^ 

From  Figure  (3.5b)  it  is  seen  that 

*  *  * 

t,  -  =  u,^  cos(})^  -  w  sin(()^  (3.51a) 

•      •  • 

n,  =  -w  cos(})^  -  u^   sincf),  (3.51b) 
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Fig.    3.5 
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similarly, 

^22  ~  ^20  *^°^^2  ~  ^  sin(()2  (3.51c) 


n^  =  -w  cos(})^  -  u_-,  sin(|)p  (3.51d) 


Using  small  angle  approximation  and  neglecting  mixed 
products  of  the  displacements  one  finds  the  components 


^11  =  ^21  =  u,o  +  ^  (3.52a) 


n.  =  -w  (3.52b) 


*11  =  *21  =  -i^  -  ^  '3-52C) 


plus  five  other  components  which  can  be  obtained  by 
permutation  of  the  indices. 

The  difference  between  the  external  and  internal  energy 
rates  must  be  zero.   If  the  procedure  which  led  to  equation 
(3.49)  is  applied  to  all  terms  of  equation  (3 . 47a)  and  the 
results  subtracted  from  the  external  work  rate  (3.47b) 
formed  with  the  aid  of  equations  (3.52)  ,  the  final  result 
must  equal  zero. 

Since  the  displacement  rates  are  arbitrary,  the  coefficient 
of  each  displacement  rate  must  be  identically  zero.   The  line 
integral  terms,  with  two  exceptions,  simply  show  that  the 
internal  forces  and  moments  at  the  boundary  are  equal  to  the 
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external  forces  and  moments  at  the  boundary.   Since  Q   =  Q 
at  the  boundary,  the  line  integral  terms  which  relate  Q   to 
the  internal  moments  at  the  boundary  are  presumed  to  be  valid 
everywhere  within  the  shell 


°1  =  k    [«2^1  -  <V,1  -  <"l2^>,2  -  ^12^,2'     '3.53a) 


Q2  =  h    I"l^,2  -  '^2^', 2  -  <«12S',1  -  "l2^1^     <3-"b): 


The  area  integral  terms  lead  to  the  following  results 


(M^B)  ^    M„B  ^     (M^^A)  ^    MA 


AB(P2  -  yii^o)  -  N^A^2  -  ^Nl2^^1  ^  ^N2^^2  "^  ^12^,1  ^ 

(MA)       M  A  _     (M,  _B)       MB 

f l£. L^iA   +  ii Li  +   -^;  ^-^  =  0  (3.53d) 

^2         2         ^2  2 


N  AB    N  AB     N  Bw  N  Aw 

AB(P3  -  ,w)  -  -1--  4--.  (-V^)a  ^  (-V^\2  ^ 

(M.B)  (MA) 

^^12-,2\l  ^  ^^12-,l\2  ^  f-V^^l  ^  [-V^^2  ^ 

^12^  2         ^12^  1  ^1^  2 

MB 

^"A^^\l  "  °  (3.53e) 
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Retention  of  the  nonlinear,  finite  deflection,  terms  in 
the  strain  rate  equations  leads  to  the  terms  involving  the 
transverse  deflection  w  in  equation  (3.53e),  which  is  the 
equilibrium  equation  in  the  Z  direction.   These  terms  do  not 
appear  in  the  case  of  infinitesimal  deflections. 

The  infinitesimal  deflection  equilibrium  equations  of 
this  thesis  have  been  compared  with  those  of  several  authors 
for  shells  of  revolution  and  have  been  found  to  be  in  exact 
agreement. 
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CHAPTER  4 
THIN  SHELL  MATERIAL  RELATIONS 

Section  1.   Yield  Surfaces 

It  was  noted  in  Chapter  2  that  the  yield  surface  and 
the  associated  flow  rule  for  rigid-perfectly  plastic  thin 
shells  are  usually  described  in  terms  of  the  generalized 
stresses  and  generalized  strain  rates  defined  in  Chapter  3. 

One  of  the  major  assumptions  of  shell  theory  is  that  the 
transverse  normal  stress  is  negligible,  and  that  the  state 
of  stress  at  a  point  on  any  parallel  surface  of  the  shell  is 
essentially  biaxial.   The  state  of  stress  can  be  expressed, 
then,  in  terms  of  two  principal,  in-surface  stresses,  a, 
and  a^.   The  yield  surface  of  a  biaxial  state  of  stress 
expressed  in  terms  of  the  principal  stresses  was  seen  to  be 
a  plane,  closed,  convex  curve  (Figure  2.1). 

Substitution  of  the  principal  stresses  into  the 
generalized  stress  definitions  (equations  3.45)  produces 
a  yield  criterion  in  terms  of  four  principal  generalized 
stresses,  two  membrane  forces  (N,  and  N„)  and  two  bending 
moments  (M   and  M„) .   The  corresponding  yield  surface  is 
described  in  four-dimensional  space,  as  is  the  generalized 
strain-rate  vector.   A  generalized  stress  yield  point  is 
reached  and  generalized  plastic  strain  can  begin  at  a  shell 
cross  section  only  when  the  yield  condition  for  stress  at  a 
point  has  been  reached  for  every  point  of  the  cross  section. 
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That  is,  plastic  flow  at  a  cross  section  is  constrained  if 
any  level  of  the  cross  section  remains  rigid. 

No  mention  has  been  made  of  the  transverse  shear  forces, 
Q.   Since  the  shell  theory  used  here  assumes  that  there  is  no 
transverse  shear  strain,  it  is  assumed  that  the  transverse 
shear  stresses  and  the  resultant  transverse  shear  forces 
have  no  effect  on  yielding.   The  transverse  shear  forces  are 
treated  as  reaction  forces  necessary  for  equilibrium. 

If  the  principal  generalized  stresses  are  not  known, 
the  yield  criterion  will  be  expressed  in  terms  of  the  six 
independent  generalized  stresses  (N  ^   and  M  „),   The  yield 

CXp  Otp 

surface  is  then  described  in  six-dimensional  generalized 
stress  space. 

The  class  of  problems  consisting  of  axisymmetric 
loading  and  response  of  thin  shells  of  revolution  is  one 
for  which  the  principal  stress  directions  correspond  to  the 
coordinate  directions.   The  four-dimensional  generalized  stress 
yield  surface  is,  therefore,  directly  applicable.   The  yield 
surface  based  upon  the  Tresca  yield  condition  was  developed 
by  Onat  and  Prager  [4-1] .   The  VonMises  yield  surface  in 
four-dimensional  generalized  stress  space  vjas  developed  by 
Hodge  [4-2].   To  the  author's  knowledge  neither  yield 
condition  has  as  yet  been  extended  to  six-dimensional 
generalized  stress  space. 

The  VonMises  yield  condition  for  a  shell  of  revolution 
is  considerably  more  complex  for  analytical  work  than  the 
Tresca  condition.   Onat  and  Prager 's  Tresca  condition  is 
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considered  in  this  thesis  to  be  the  exact  yield  condition. 

For  problems  in  six-dimensional  stress  space, and  often 
in  four-dimensional  stress  space  problems  as  well,  approximate 
yield  surfaces  are  used.   These  simplified  yield  surfaces  may 
be  necessary  in  order  to  make  the  resulting  equations  tract- 
able.  Although  the  results  obtained  from  approximate  yield 
conditions  will,  in  general,  be  in  error,  the  bounding  theorems 
of  limit  analysis  allow  the  exact  results  to  be  bounded  with 
two  approximate  yield  surfaces  of  the  same  shape,  one  of  which 

circumscribes  the  exact  surface  and  one  of  v/hich  inscribes 
the  exact  surface  [4-3] . 

The  two-moment  limited  interaction  yield  surface  proposed 
by  Hodge  [4-4]  and  shown  in  Figure  4.1,  for  example,  circum- 
scribes the  Tresca  yield  surface,  while  a  similar  surface 
with  a  scale  0.618  times  that  of  the  original  inscribes  the 
Tresca  surface.   The  results  which  would  be  obtained  through 
use  of  the  Tresca  surface  must  lie  between  the  results  obtained 
through  use  of  the  two  limited  interaction  surfaces. 

A  two-direction  limited  interaction  yield  surface 
suggested  by  the  Tresca  yield  surface  of  a  curved  beam  [4-5] 
and  shown  in  Figure  4.2  also  circumscribes  the  exact  Tresca 
surface.   The  inscribing  surface  also  has  a  scale  0.618  times 
that  of  the  circumscribing  surface  (Appendix  A) . 

Closer  bounding  may  be  possible  if  the  stress  points  of  a 
given  problem  lie  only  on  a  small  portion  of  the  yield  surface. 
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Section  2.   Discontinuities 

In  the  mathematical  theory  of  rigid-perf ectly  plastic 
shells,  the  possibility  of  discontinuities  in  displacement, 
displacement  rate,  and  stress  fields  exists.   It  should  be 
recognized  of  course  that  what  may  be  a  discontinuity  in  the 
rigid-perf ectly  plastic,  midsurface  representation  of  a  shell 
would  in  the  actual  shell  be  a  narrow  zone  across  which  the 
"discontinuous"  field  quantity  varies  rapidly. 

Mathematical  discontinuities  and  the  associated  notion 
of  a  generalized  hinge  are  useful  ideas  in  the  model  of  an 
ideally  plastic  shell.   Not  all  discontinuities  are  acceptable, 
hov/ever,  and  those  which  violate  the  assumptions  of  the  shell 
theory  should  be  identified  and  eliminated. 

This  examination  extends  the  work  of  Hopkins  on  flat 
plates  [4-7]  to  shells. 

First  let  us  examine  displacements  and  displacement 
rates. 

The  in-surface  displacements  and  displacement  rates 

(u  ^  and  u  ^,  a=l,2)  may  be  discontinuous.   Fracture  is  not 
aO       aO      '    '         J. 

considered  in  this  thesis,  and  these  discontinuities  are 
interpreted  as  severe  local  deformation. 

A  discontinuity  in  the  transverse  displacement  or 
displacement  rate  (w  or  w)  would  correspond  to  severe  trans- 
verse strain.   The  basic  Love-Kirchof f  assumption  requires 
transverse  strain  to  be  zero.   Therefore,  w  and  w  are  assumed 
to  be  continuous. 
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The  continuity  of  a  quantity  across  a  line  in  a  shell 
imposes  certain  requirements  on  the  derivatives  of  that 
quantity. 

Figure  4.3a  shows  a  simple  curve  F (t)  in  the  midsurface 
of  a  shell  element.   The  s  and  n  directions  are  tangent  and 
normal  to  the  curve,  respectively.   The  curve  T    is  shown 
in  Figure  4.3b  passing  through  a  differential  element  of 
length  ds  and  width  dn  centered  on  the  point  P(t).  T   divides 
the  shell  element  into  two  regions  as  shown,  R   and  R  .   The 
points  P',  P',  P^,  P",  P",  P",  P    ,    and  P'  lie  on  the  boundary 
of  the  differential  element. 

Let  us  define  the  discontinuity  of  a  field  quantity  G 
as  one  crosses  the  curve  T    at  the  point  P  as 


[G]p  =  G^  -  G  (4.1) 


where 


G   =  lim  Gp„ 


pp"  ->■  0 

and 


G   =  lim  G  , 

pp'  ^  0 
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Fig.    4.3  (b) 
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ra^di) 


(c) 


Fig.    4.3     (con't) 
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Let  us  suppose  that  [G]  =  0  all  along  the  curve  T,   Then 


[G]^    =  [G]p  =  [G]p   =  0  (4.2) 

^1        ^        2 


Therefore, 


lim     [G]p   -  [G]p 


That  is,  there  is  no  discontinuity  across  T   of  the  rate  of 
change  of  G  along  T. 

If  we  look  at  the  position  of  F  at  two  instants  of  time, 
t  and  t+dt,  (Figure  4.3c) 


f<^ip(t)  =  "5ip(t+dt)  =  °  '*■*> 


Since  the  curve  T    is  assumed  to  be  moving,  the  change  in 
the  discontinuity  as  a  function  of  time  is  a  function  of  both 
position  and  time,  and  leads  to  a  substantial  or  material 
derivative  [4-6] . 


lin,    [Glp(     )  -  tG]p(^  ^    ,         ^  [|?]  =  0        (4.5) 
dt-0         "^t  '"     ^^ 


where  V  is  the  instantaneous  velocity  of  T   normal  to  the 
curve. 
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Equation  (4.5)  presents  the  following  three  possibilities: 


If  V  =  0 


(|f]  =  0 

[^5—]  is  not  necessarily  zero  (4.6a) 


If  V  ^  0 

either 

or 


[|f]  =  [|f]  =  0  (4.6b) 


[|f]  =  -i[|f]  (4.6C) 


Applying  equations  (4.3  and  4.6)  to  the  continuous 

3w, 
dt' 


•  •    3w» 

quantities  w  and  w  one  finds  (recalling  that  w  =  -r^) 


[||]  =  0  (4.7a) 


[|^]  =  0  (4.7b) 

dS 


If  V  =  0   (stationary  hinge) 


[■r— ]  not  necessarily  zero  (4.7c) 

dn 


[t: — ]  not  necessarily  zero  (4.7d) 

dn 


[w]  =  0  (4.7e) 

70 


If  V  7^  0    (travelling  hinge) 


[|^]  =  0  (4.7f) 

an 


and  either 


or 


[|^]  =  [w]  =  0  (4.7g) 

on 


tf^l  =  -itS)  (4.7h) 


The  requirements  for  continuity  of  forces  and  moments 
across  the  line  V    at  the  point  P  are  obtained  by  looking 
at  the  equilibrium  of  the  differential  element  of  Figure  (4.3b) 
This  element  is  redrawn  in  Figure  (4.4)  with  the  membrane  and 
transverse  shear  forces  indicated. 

The  condition  of  vertical  force  equilibrium  can  be  ex- 
pressed as 


p  It         p "  p '         p " 

1  2  2  2 


(    /      -    /    )    Q      cosco    dn   +    (    /      -    /    )    0      cos5    ds 
pi         pi  s  pi         p  II         n 

1  2  1  1 


p  I  p  n  p "  p  II 

2  2  1  2 

(    /  +   /    )  N      sin6    ds    -    (    /  +    /    )    N      sinco      dn 

p  I  p"  ^                                 p  I  p  I         ^ 

1  1  ^1  ^2 
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Qs     ^f^ 


Nc 


Fig.    4.4 
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/   /      (P-3  "  VJw)  dnds  =  0  (4.8a) 

pi  p"  p"  p«    ■^ 
^1   1   2   2 


If  we  now  allow  P'  and  P"  to  approach  P^  and  P'  and  P" 
to  approach  P^,  dn  and  d(})  approach  zero.   Then 

^2 

I   [Q„]  ds  =  0  (4.8b) 

Finally,  allow  P   and  P   to  approach  P  to  find  that 


[Qj^lp  =  0  (4.8c) 


Similar  operations  with  the  force  and  moment  equilibrium 
requirements  in  the  s  and  n  directions  yield  the  following 
results. 


[M^]^  =  m^J^   =  [N  ]p  =  [N  ,]p  =  0  (4.8d) 

n  P      ns  P     n  P     ns  P 


The  conclusion,  therefore,  is  that  only  Q  ,  N   and  M 
may  be  discontinuous  across  the  curve  F. 
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CHAPTER  5 
RESPONSE  OF  ARBITRARY  THIN  SHELL  TO  DYNAMIC  LOAD 

Experience  has  shown  that  bodies  which  have  been  plas- 
tically deformed  often  have  regions  with  little  or  no 
permanent  strain  separated  by  relatively  narrow  regions  which 
have  been  severely  strained.   Although  the  regions  of  severe 
deformation  have  a  finite  width,  they  are  often  narrow  enough 
to  suggest  the  appearance  of  a  line.   These  lines,  or  hinges, 
are  normally  thought  of  in  conjunction  with  severe  bending, 
which  causes  a  discontinuity  in  the  slope  of  the  body.   It 
was  noted  in  Chapter  4,  however,  that  discontinuities  in 
other  quantities  can  also  occur  in  a  body.   The  term  "hinge" 
in  this  thesis  will  be  construed  to  mean  any  line  across 
which  an  allowable  discontinuity  in  the  stress  field  or  the 
displacement  field  (or  their  derivatives)  exists. 

Although  line  hinges,  that  is,  hinges  of  zero  width, 
do  not  actually  occur,  the  representation  of  such  a  hinge 
can  be  mathematically  consistent  within  the  framework  of 
the  thin  shell  assumptions  so  far  introduced. 

Since  hinges  appear  to  be  lines,  one  would  expect  the 
dissipation  of  energy  at  a  line  hinge  to  be  represented  by 
a  line  integral.   Such  was  the  representation  of  the  external 
work  rate  at  a  shell  boundary  used  in  Chapter  3  (equation 
3.47b).   It  should  therefore  be  instructive  to  look  for  line 
integral  representations  of  shell  energy  dissipation  expressions 

If  there  is  no  dissipation  of  energy  at  the  shell  boundary 
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(a  restriction  which  will  not  be  maintained  for  long) ,  the 
external  energy  dissipation  rate  for  a  shell  is  the  work,  rate 
of  the  external  forces  acting  on  the  shell. 


Dg  =  /HP^  -  yu^o^^lO  ■"  ^^2  -  ^^20^^20 


+  (P^  -  yw)w}dS  (5.1) 


Each  expression  within  parenthesis  may  be  replaced  by  an 
expression  involving  the  generalized  stresses.   The  resulting 
equation,  which  will  not  be  reproduced  here,  is  simply 
obtained  by  use  of  equilibrium  equations  (3.53c,  d  and  e) . 

Each  term  of  the  resulting  expression  can  be  integrated 
by  parts  over  any  region  within  which  the  term  is  continuous. 
This  process  leads  to  a  line  integral  along  the  curve  bounding 
the  region  and  an  area  integral  within  the  region. 

Let  us  carry  out  the  integration  of  a  typical  term,  the 
first  term  in  the  u,-  integrand,  over  a  region  of  the  mid- 
surface  bounded  by  the  curve  C  (Figure  5.1). 


/[/N2B^^u^QdnJdn2  =  /{N2B  u^q}  dn2  -  ffBil^^u^^)  ^^dr^^dr,^ 


^2^1  ^1 


(5.2) 


where  {}    indicates  that  the  bracketed  quantity  is  to  be 

^1 
evaluated  at  the  n-.  limits.   For  a  given  curve  C  the  values 

of  n-,  and  r\      are  not  independent,  and  the  n-i  limits  must  be 
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Fig.     5.1 
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expressed  in  terms  of  n^  in  order  to  continue  integration  of 
the  first  term  on  the  right  hand  side  of  equation  (5.2). 
This  integration  process  is  indicated  schematically  in 
Figure  (5.1).   The  narrow  strip  indicates  the  integration 
over  ri-,  and  the  arrows  indicate  the  r\y   integration.   The  t]^ 
integration  can  be  thought  of  as 


/{N2Bu^Q}dn2    -    /{N2Bu^Q}dTi2 
(Path    RQS)  (Path   RPS) 


The  sign  of  the  second  term  is  changed  if  the  integration 

path  is  reversed  to  SPR.   The  entire  integration  is  then  in 
the  counter-clockwise  direction  around  the  closed  curve  C 
and  may  be  represented  as  a  line  integral.   A  counter- 
clockwise integration  path  is  considered  positive  by 
convention  [5-1] .   The  same  reasoning  applied  to  an  integral 
of  the  form 


/{F)^^  dn, 


leads  to  a  clockwise  integration  path  which,  by  convention, 
is  considered  negative. 

The  final  result  for  the  dissipation  term  being  considered 


is 


//N2B^^u^Qdn^dn2  =  ^N2Bu^Qdn2  -  //B(N2U^Q)^^dn;L^n2   (5.3) 
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The  closed  curve  C  would  be  chosen  as  a  line  across 
which  a  discontinuity  in  (N„Bu.q)  exists.   There  may  be 
several  curves  of  discontinuity  for  each  term  of  the 
dissipation  integral.   The  contribution  from  each  term  is 
the  sum  of  its  associated  area  and  line  integrals  from  (and 
around)  each  region.   Not  all  terms  would,  in  general,  have 
discontinuities  at  the  same  curves.   Each  term  considered 
might  have  its  own  particular  regions  of  continuity  and 
hinge  lines.   Therefore  each  term  should  be  separately 
summed.   For  simplicity  in  presenting  the  results  of  converting 
the  dissipation  integral  to  the  sum  of  line  and  area  integrals, 
only  one  summation  index  will  be  used,  merely  to  serve  as  a 
reminder  that  summation  of  the  contribution  from  the  different 
regions  is  required.   If  the  procedure  described  above  is 
applied  to  each  term  of  the  dissipation  integral  which  can 
have  an  allowable  discontinuity,  the  final  result  is 


MM  M 

n  11  2 


N^Bw'  T  (M^B)  T  2Mt^A      MB     . 

^   A        ^12  ,2       A       ''l2,2       A        ^1 


+  (N^w  +  M^)|w^^  +  (N^2^  +  Ml2^^',2^^^2 
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MM  M 

(5{[N^   -   "2   -   r|  ^   R^l*^20    -    2<Ni2   ^   nf  )^^10 


N„AW 


(M^A) 


2M, ^B    ,  M,A 


r_J l1   +    n      w         + '—  +   M  +    — — —   -       •*-     '^]    w 

^       B  12     ,1  B  12,1  B  B       ^ 


+     (N^w    -f   M^)^   w^2    -^     ^N^2^   -^   ^2^    ^,1^^"^! 


^10.  _^2^10, 


^//{-bcn^u^q)^^  +  N^Bu^oa  ^  ^l^(R^^l  -  ^(^^^l 


u 


20, 


■^   ^12^^20,1    ^   ^^^12^20^1    -^   ""ll^^-R^Kl 


^12''20 
^   ^(-^^l    ^   ^12^^10,2    ^   ^^^12^10^2 


"^10  ^12^^10 

^   ^12^^Rf\2    ^   ^^^l^^2    -   ^(V20^2    ^   ^2^^20,2 


u  M   u  N  N 

20  120  1  2* 

^   «2^'^',2   ^  ^<-^',2   ^  ^^%  *    rJ)"   ^   »12,2«,1 


2M, _A    _W  M_B    ,W    .  2M, ^B    ,W    _ 

+   M  w         +         12    ,2    ,1    _      2    ,1    ,1  12    ,1    ,2 

''l2,l    ,2  A  A  B 


M   A      W  N    BW  W 

B       '       -    '^A^>,1  "   -'"12",2',1  "    -   V<^',1 


N   AW  W 
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If  the  external  energy  dissipation  rate  expressed  in 
terms  of  external  load  and  displacement  rates  (equation  5.1) 
were  equated  to  the  external  dissipation  rate  expressed  in 
terms  of  the  generalized  stresses  and  the  displacement  rates 
(equation  5.4),  and  if  the  generalized  stresses  were  expressed 
in  terms  of  the  displacements  and  the  displacement  rates 
through  use  of  the  material  yield  surface  and  the  rule  of 
normality,  then  the  displacement  field  could  be  found  as  a 
function  of  time,  in  principle.   The  extreme  complexity  of 
such  an  endeavor  has  precluded  exact  solution  of  such  a 
problem  except  for  very  simple  cases. 

In  the  next  chapter  the  general  results  of  this 
chapter  will  be  specialized  to  a  shell  of  revolutiorv  and  an 
approximate  solution  of  the  shell  response  to  dynamic  load 
will  be  presented. 

It  should  be  noted  that  the  general  dissipation  equation 
(5.4)  can  include  hinge  lines  which  coincide  with  the  physical 
boundaries  of  the  shell.   Therefore,  the  external  work  done 
at  a  shell  boundary  can  also  be  accounted  for,  effectively 
removing  the  restriction  imposed  before  equation  (5.1). 
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CHAPTER  6 
AN  APPROXIMATE  METHOD  FOR  DETERMINING  THE  RESPONSE 
OF  A  THIN  SHELL  OF  REVOLUTION  TO  AXISYMIIETRIC  DYTJAMIC  LOAD 

In  Chapter  5  the  external  energy  dissipation  rate  of  a 
general,  thin  shell  of  uniform  thickness  was  developed  in 
terms  of  area  and  line  integrals  involving  the  generalized 
stresses,  the  displacement  rates,  the  transverse  displace- 
ment (if  deflections  are  finite),  and  their  derivatives. 
The  line  integrals  represent  mathematically  the  possibility 
of  hinges  in  a  rigid-perfectly  plastic  material. 

In  this  chapter  the  general  results  indicated  in 
Chapter  5  are  specialized  to  the  case  of  a  thin  shell  of 
revolution  which  is  subjected  to  an  axisymmetric  dynamic 
load.   Neither  buckling  nor  nonsymmetric  response  is 
considered. 

A  commonly  used  coordinate  system  for  shells  of 
revolution  is  indicated  in  Figure  6.1  along  with  the  positive 
directions  of  transverse  and  meridional  load  and  displacement. 
It  is  assumed  that  there  are  no  loads  or  displacements  in  the 
circumferential  direction. 

Table  6.1  relates  the  symbols  used  for  general  shells 
in  previous  chapters  to  the  symbols  commonly  used  for  this 
axisymmetric  class  of  problems. 
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Fig.    6.1 
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<J> 

e 

^1 

^2 

\ 

R^sin4) 

TABLE  6.1 

AXI SYMMETRIC  SHELL  PROBLEM  SYMBOLS 
AND  THEIR  GENERAL  SHELL  PROBLEM  COUNTERPARTS 


General 
Z 

n 

2 

1 
'2 


'1 

n, 
^] 

R, 

A 

B 

^10 
w  w 

^(|)  ^10 

^4)  ^1 

^<t>e  ^12 

\'  %'•  Q^  ^1'  ^1'  Qi 

Mq;    Nq;    Qq  m^;    n^;    Q2 

p  p 

<}>  1 

p  p 

3  3 
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since  the  geometry  of  this  class  of  shells  does  not 
vary  with  6,  the  derivatives  of  all  geometric  quantities 
with  respect  to  6  are  identically  equal  to  zero. 

Because  of  the  assumed  axisymmetry  of  the  load  (and  the 
shell  response)  the  derivatives  of  all  stresses  and  dis- 
placements with  respect  to  9  are  also  identically  equal  to 
zero. 

The  non-zero  rate  equations  (3.44a,  b,  d,  e)  reduce  to 


£.  =  ^  (u'  +  w  +  V^)  <6-l^) 

Ir,   =  h      (ucotd)  +  w)  (6.1b) 

e   R2 


1    .     .     (w'-u)R' 

K  .  =  1^  [u'  -  w"  +  -— ^]  (6.1c) 

(p         R^  R^ 


cot*   / •   •  t  »  /  r      1  j\ 

Kg  =  p-^  (u-w')  (6. Id) 


where 


The  equilibrium  equations  (3.53)  reduce  to 


R^R^  sind)(P^  -  yii)  -  N^R^  cos*  +  (N^R^  sin*)'  + 
12         (p  6  1  w    Z 


(MR   sin(}))  ' 
— ^^ -  M-cos(})  =  0  (6.2a) 


\ 
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Rt  R„    sin(t)(P-    -    yw)    -   N.R^    sine})   -   N^R,    sincj)   + 
J.    z  J  9    i^  d    J. 


N   R      sin(})  w'                 (^A^o    sin(}))  ' 
i^^ )•    +    [      '^    \ ]•    -    (MqCOSc)))'    =    0       (6.2b) 


M' 


Noting  that  the  differential  midsurface  area,  dS ,  is 


dS  =  (27tR  sine}))  R^d(}) 


and  neglecting  external  work  at  the  shell  boundary,  the 

external  energy  dissipation  rate  can  be  expressed  in  terms 
of  external  loads  as 


Dg   =    2tt/[(P      -    yu)u   +    (P^    -    yw)w]    R^R^    sin(j)dc!)  (6.3) 

The  dissipation  rate  in  terms  of  area  integrals  and 
line  integrals  may  be  obtained  from  the  general  results 
(equation  5.4)  or  by  applying  the  procedure  of  Chapter  5  to 
equation  (6.3).   Application  of  the  procedure  to  a  portion 
of  the  dissipation  integral  is  demonstrated  in  this  chapter. 

Symmetry  limits  hinge  curves  in  the  rigid  perfectly- 
plastic  material  to  circles  (cj)  =  constant),  and  the  regions 
of  continuity  for  each  term  of  the  dissipation  integral  are 
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therefore  each  bounded  by  two  circles.   Since  there  is  no 
variation  with  6,  the  problem  is  essentially  one-dimensional. 

Regions  in  which  a  typical  term  of  the  dissipation 
integrand  (expressed  in  terms  of  the  generalized  stresses) 
is  continuous  are  indicated  in  Figure  6.2. 

Integration  by  parts  of  such  a  typical  term  gives 

/(-N.R_sin(}))  •  u  d(|)  =  E/N^R^sin(})u'dd)  -  E{N^R_sin(J)u}  (6.4) 
4)2  ^^(j)2        ^         ^      (^    2 

where  {  }  indicates  that  the  expression  is  to  be  evaluated 

between  the  upper  and  lower  limits  ^      corresponding  to  the 

nth  region. 

Indicating  the  value  of  (N.R^  sin(j)u)  evaluated  by 

approaching  6   from  below  as  (  )   and  from  above  as  (  )  , 
^'-  ^  ^m  m   .  m 

the  second  term  on  the  right  hand  side  of  equation  (6.4)  can 
be  expressed  as 


-  Z  {N.R_sin(|)u}   =  -  E  [(  )     .^    "  (  )^]  (6.5) 

<n    2.  m+1       m 

n    ^  m 


If  the  summation  over  m  is  expanded  and  the  terms 
grouped  according  to  the  pLastic  hinge  index  m,  this 
becomes 


E  {N  R^  sin({)u}  =  -[(  )"  -  (  )^]  -  [(  )2  -  (  )\'^    " 
n    ^ 


+  E  [N^R^sin4)u]  (6.6) 

d)  2       m 
m    ^ 

86 


Axis   of 


Symmetry 


fr]~l 
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<»■ 


where  [F]   indicates  the  discontinuity  in  F  across  the  hinge 

circle  (j)  =  6  .   That  is,  [F]_  =  F^  -  F~. 
m  m    m    m 

If  attention  is  restricted  to  shells  with  smooth 

generators,  geometrical  factors  such  as  R^sincJ)  are  not 

discontinuous.   Let  the  notation  (G)   indicate  that  the 

m 

function  G  is  to  be  evaluated  at  d)  .   Then  equation  (6.4) 

m 

can   be   written    as 


/(-N   R      sin(J))  •    ud(j)    =      ^/N,R2    sintf)   u'    dc|) 
<P  n      ^ 


+    S(R„    sin(|))^    [N.u]^  (6.7) 

2  mom 

m  ^ 


Investigation  of  discontinuities  across  the  hinge 
circles  based  on  the  assumption  that  the  w  and  w  are 
continuous  everywhere  and  based  on  the  interpretation  that  a 
mathematical  discontinuity  in  u  merely  represents  a  very 
large  local  strain  yields  the  following  results  (c.f. 
Chapter  4) . 


[M,]  =  [N.]  =  [Q.]  =0  (6.8a) 

(p  (^  (p 


[w']  may  be  different  from  zero  at  a 

stationary  hinge.  (6.8b) 

[w']  may  be  different  from  zero  at  a 

travelling  hinge.  (6.8c) 
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since  [M^]  and  [Q,]  are  both  zero,  equation  (6.2c)  shows 

that 

R^  cot(J) 


With  the  discontinuity  results  of  equations  (6.8)  and 
(6.9)  and  the  notation  of  equation  (6.7),  the  external 
dissipation  rate  expression  for  the  axisymmetric  case  may 
be  written  as 


D  -M  ^'^^^2   sine}). 

1  9   1  ml 


M,R„  sin(}) 
+  /.  {ZN,R_  sine})  -  ZN^R_  sin(|)  +  E-  ''^ 


.<.    2   ^^'v    r  e  2    ^         R 
^  n  ^  1  m 


1 


MR      sin(})    . 
^-| }u'dcf)    +/x^N    R^    +   NgR^}sin4)   wd<j) 


R'R      sin(}) 

^k1  ^  ^K^  1 


R^    sine}) 

-    M„cos(^}w'd(|)    -/^{Z-^^ (N^w   +   M.)}w"d(|) 

e  *    K         1  ^  ^ 


M. 
+    ^{R^sinMN      +Rf)>q[u]       "    MR2Sin<J))        [N^u] 
q  l^^p  ^  ^ 


R^sincJ)  R   sin(}){N   w   +  M    } 

y       Ri      y       e     y       ^  Ri  K  K 
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+  K-^ ^)^  [W]^  (6.10) 


It  should  be  noted  that  the  nonlinear  terms  in  w  appear 
only  if  finite  deflections  are  considered.   One  of  these,  the 
final  term  of  equation  (6.10),  can  be  shown  to  be  zero. 

Continuity  considerations  revealed  that  [w']  could  be 
different  from  zero  at  a  stationary  hinge  only.   Now  [w']  ^  0 
implies  that  w"  and  w"  are  infinite  and,  therefore,  that  k. 
is  infinite.   For  the  exact  Tresca  yield  surface  of  shells 
of  revolution  (Onat  and  Prager  [6-1]),  normality  requires 
the  ratio  K^/e,  to  be  finite  unless  N,  is  zero.   Finite 
deflections  give  rise  to  membrane  forces  and,  therefore, 
e,  would  also  be  infinite.   If  the  shell  stretched  at  a 
hinge  with  an  infinite  rate  for  a  finite  period  of  time, 
infinite  thinning  would  occur.   The  finite  deflection  teann 
involving  [w']  is  therefore  zero. 

Solution  of  equation  (6.10)  in  terms  of  displacements 
as  functions  of  time  and  position  is  a  very  complex  problem, 
and  is  made  more  so  because  of  the  presence  of  travelling 
hinges. 

In  order  to  simplify  the  problem,  let  us  restrict  our 
attention  to  problems  with  transverse  loading  only. 

For  the  remainder  of  this  chapter  it  will  be  assumed  that. 
P.  =  0  and  that  ii  =  0.   Therefore,  the  integrals  and  hinge 
terms  involving  u  are  identically  zero. 
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Another  approximation,  suggested  by  Sawczuk  [6-2]  and  by 
Jones  [6-3]  for  the  static  and  dynamic  response,  respectively, 
of  flat  plates,  is  also  used. 

Although  the  discontinuity  [w']  can  only  exist  at  a 
travelling  hinge,  let  us  treat  such  a  hinge  as  if  it  is  in 
fact  stationary  and  formed  at  the  same  location  at  which  a 
static  collapse  hinge  would  form.  Identical  displacement 
profiles  for  static  collapse  and  for  dynamic  response  have 
been  found  for  cases  such  as  square  plates  [6-4] ,  circular 
plates  [6-5]  ,  and  cylindrical  shells  [6-6]  ,  if  the  applied 
load  does  not  greatly  exceed  the  static  collapse  load. 

If  the  external  dissipation  rate  expression  for  the 
transverse  load  problem  in  terms  of  generalized  stresses  is 
equated  to  the  dissipation  rate  expression  in  terms  of 
external  load,  the  result  is 


^r-  ^•^A'fN.R-    sin<t>   +  N.Rt    sin(})}wd4)    + 

277992  61  -r-r 


lU^  ^"\r"    *V  ^  V  '  '^''2  ^i"*''  h^  "''^* 


/,MqCOS(J)   w'dd)    -    If,    2  /-M         ,    r,    \     *njA 

d)    6         ^  ,     4) :?; (N,w   +   M.)    w   dci 

k    ^         R,  (})  (f)  ^ 

Rp    sint}) 

Z{ — 5 (N.w   +   M , )  }     [w']       =    /.R.R^    sinct)  (P^-pw)wd(})     (6.11) 

R.,  9  9K  K  9I2  3 

where  boundary  loading  has  been  omitted  for  simplicity. 

91 


To  use  equation  (6.11)  in  estimating  the  permanent 
deflection  of  a  shell  of  revolution  subjected  to  an 
axisymmetric,  transverse  dynamic  load,  one  begins  by 
assuming  a  kinematically  admissible  velocity  profile 
whose  shape  is  time  independent.   The  static  collapse  profile 
is  suggested,  if  known. 


w{(}),t)  =  w^Ct)  -f  ((}))  (6.12) 


where  f(cj))  is  the  deformation  velocity  mode  and  w  (t)  is 
the  velocity  mode  amplitude. 

The  velocity  profile  of  the  assumed  mode  of  deformation 
is  used  with  equations  (6.1)  to  determine  the  strain  and 
curvature  rates. 

The  rate  equations  are  used  in  conjunction  with  the 
assumed  yield  surface  and  its  associated  flow  rule 
(normality)  in  order  to  determine  the  values  of  the 
corresponding  generalized  stresses. 

The  expressions  derived  for  the  generalized  stresses 
are  entered  into  equation  (6.11)  and  the  indicated 
summations  and  integrations  performed.   The  result,  after 
factoring  out  common  coefficients,  is  a  second  order 
differential  equation  in  w   and  time.   The  equation  will  be 
nonlinear  if  finite  deflections  have  been  retained. 
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The  exact  formulation  of  this  differential  equation  is 
dependent  upon  several  factors  such  as  shell  geometry/  assumed 
deformation  mode,  and  assumed  yield  surface.   Specific 
examples  are  given  in  Chapter  8. 

The  method  of  solving  the  differential  equation  depends 
upon  its  formulation  and  upon  whether  or  not  the  pressure 
pulse  shape  is  known.   If  the  shape  of  the  pulse  is  not 
known  or  if  it  is  a  high  pressure,  short  duration  pulse,  the 
problem  can  be  treated  as  impulsive  loading. 

If  the  time  history  of  the  pressure  pulse  is  to  be  used, 
the  initial  conditions  are 


WQ(t=0)  =  WQ(t=0)  =  0 


If,  however,  impulsive  loading  is  assumed,  the  initial 
conditions  are 


WQ{t=0)  =  0 


WQ(t=0)  =  vVq 


where  V   is  the  actual  impulsive  velocity  at  the  location  of 
w  ,  and  V  is  the  Martin  and  Symonds  mode  factor  for  optimum 
matching  of  the  actual  impulsive  problem  with  the  mode 
approximation  (time  independent  profile)  [6-7].   Equation  (6.11) 
reduces  to  the  flat  plate  approximation  of  Jones  [6-3]  . 
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CHAPTER  7 
APPROXIMATE  RESPONSE  OF  NON-SYMMETRIC  SHELL  TO  DYNAMIC  LOAD 

The  solution  technique  of  Chapter  5  can  be  applied,  in 
principle,  to  any  shell  which  satisfies  the  usual  thin  shell 
theory  assumptions.   In  this  chapter,  however,  an  approximate 
technique  is  investigated  which  is  not  derived  from  the 
previous  results.   The  technique  for  dynamic  loads  developed 
here  is  a  simple  extension  of  a  static  collapse  method  pro- 
posed by  Janas  for  the  limit  analysis  of  non-symmetric  shells 
[7-11] .   Finite  deflections  have  not  been  included. 

The  basic  assumption  used  in  this  approach  is  that  the 
body  deforms  as  several  rigid  regions  separated  by  narrow 
plastic  hinges.   The  plastic  hinges  are  viewed  mathematically 
as  curves.   The  results  of  severe  membrane  strain,  twisting, 
and  bending  rates  are  viewed  in  terms  of  discontinuities  in 
velocities  and  relative  rotation  rates  between  rigid  regions 
at  the  hinge  curves.   All  internal  energy  dissipation  occurs 
in  the  plastic  hinges,  and,  according  to  the  principle  of 
virtual  velocities,  the  external  work  rate  must  equal  the 
internal  energy  dissipation  rate.   The  dissipation  rate  can 
be  expressed  as  a  function  of  the  hinge  discontinuities. 

Consider  a  hinge  separating  rigid  regions  I  and  II  to 
have  a  small  width  X,  as  shown  in  Figure  7.1.   Because  trans- 
verse shear  strain  is  assumed  to  be  zero,  there  can  be  no 
discontinuity  in  the  transverse  displacement  rate  w.   The 
only  displacement  rate  discontinuities 
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then  are  [u,^]  and  tu__]  corresponding  to  the  shell  coordinate 
axes  ri-i  and  ri^. 

The  velocity  discontinuity  ccnponents  normal  to  and 
tangent  to  the  hinge  curve  C  are  seen  to  be 


and 


[U^]  =  [u2q]  cosB   +  [u^q]  sin3  (7.1a) 


[U^]  =  [u^q]  cos3   -  [u2q]  sin3  (7.1b) 


where  3  is  the  angle  between  the  normal  to  the  hinge  curve 
and  the  coordinate  axis  T]^. 

Regions  I  and  II  are  rigid  anv^  the  length  of  the  hinge 
curve  can  not  change.   Therefore, 


e^  =  0  (7.2a) 


If  one  assumes  that  the  displacement  rates  vary  linearly 
across  the  width. 


e_  = 


[U  ] 
n 


'n     X 
and 


(7.2b) 


;  ,  =  -^  (7.2c) 


nt     X 
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By  transformation  of  axes  it  can  be  shown  that 


e-j^Q  =  — ^ —  sin3  (7.3a) 

'"20' 

^20  "^  — X cos3  (7.3b) 

[u   ]  cosB   +  [^9n^  ^^^^ 
Yq  =  -^ y ("7.3c) 


By  similar  reasoning  one  finds  the  curvature  rates  across 
the  plastic  hinge. 


^1 
K,  =  T—  sinB  (7.4a) 

^2 
^2   ~   ~~\      ^°2^  (7.4b) 

^,  cos3   +  Q^   sin3 
2k^2  =  A (^-^^^ 


where  9,      is  the  n   component  of  the  relative  rotation  rate 
between  rigid  regions  I  and  II. 

The  internal  dissipation  rate  is 


D.=/  (N£         +Ne  +N         Y+MK       + 

i        g  ^1    ^10  2    ^20  12    '0  1      1 


hinge 


"2    ^2    ■"    2«i2    '^12'    "^  "-5) 
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If  it  is  assumed  that  the  generalized  stress  values  are 
uniform  across  the  hinge  zone,  and  if  it  is  noted  that 

dS  =  AdC 

then  equation  (7.5)  can  be  written  in  the  form 


^12^^20-'  ^^^^   "*"  ^1  ^1  sin3   +  M^  ^2  cos3   + 


2M^2  ^1  cos3   +  2M^2  ^2  ^^^^^  ^^  (7.6) 


A  limited  interaction  yield  surface  is  assumed  which 
reduces  to  the  two  moment  limited  interaction  yield  surface 
(c.f.  Chapter  4)  for  N,-  =  M-^  =  0-   It  is  assumed  that 
membrane  shear  strain  can  occur  only  if  N-_  =  7  N   and  that 
twisting  can  occur  only  if  M,  _  =  7  M^,.   The  factor  of  one 
half  arises  from  the  fact  that  the  yield  stress  for  shear 
is  one  half  the  yield  stress  in  simple  tension. 

The  limited  interaction  yield  surface  is  shown  in 
Figure  7.2. 

The  scale  of  the  similar  yield  surface  which  would 
inscribe  the  exact  six-dimensional  Tresca  yield  surface  is  not 
known  because  the  six-dimensional  Tresca  yield  surface  has  not 
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been  developed.   The  scale  of  the  inscribing  surface  can  be 
no  larger  than  0.618,  which  is  the  scale  of  the  tvo  moment 
limited  interaction  surface  inscribed  in  the  Tresca  surface 
in  four-dimensional  stress  space.   The  value  0.618  is  used  in 
this  thesis  as  the  best  available  scale  factor. 

The  solution  to  a  particular  problem  would  be  carried 
out  as  follows. 

First  a  collapse  mode  is  postulated  which  involves  rigid 
body  motion  of  one  or  more  regions. 

The  rigid  body  equations  of  motion  for  each  region  are 
examined  at  each  hinge  curve  to  determine  the  velocity  dis- 
continuities and  the  relative  rotation  rates  between  regions. 

« 

The  condition  of  continuity  of  transverse  velocity  w  may  be 
used  to  determine  the  equation  of  one  or  more  hinoo  curves. 
All  velocities  may  be  expressed  in  terms  of  a  characteristic 
transverse  velocity  w^. 

The  strain  rates  and  curvature  rates  are  used  to  determine 
the  corresponding  generalized  stresses. 

The  internal  dissipation  rate  (equation  7.6)  is  equated 
to  the  external  dissipation  rate  (equation  5.1).   The  resultant 
second  order  differential  equation  in  w.  and  time  can  then  be 
solved  for  the  permanent  deformation  w. . 

The  procedure  indicated  in  this  chapter  has  been  used  to 
solve  the  problem  of  a  cylindrical  panel  loaded  impulsively 
over  a  portion  of  the  shell  inner  surface.   The  rt>:;ults  are 
reported  in  the  next  chapter. 
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Extension  of  this  technique  to  include  finite  deflections 
is  not  conceptually  difficult.   The  major  difference  is  that 
the  path  of  a  point  in  the  midsurface  of  a  shell  region  under- 
going rigid  body  rotation  must  be  recognized  as  the  arc  of  a 
circle  rather  than  assumed  to  be  the  tangent  to  the  arc. 

Rigid  body  motion  cannot  be  limited,  in  general,  such 
that  transverse  displacements  are  much  larger  than  in-surface 
displacements.   Non-linear  terms  involving  all  of  the  dis- 
placements could  therefore  be  expected  to  arise  in  finite 
deflection  problems. 
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CHAPTER  8 
APPLICATION  OF  GENERAL  METHODS  TO  PARTICULAR  CASES 

Section  1:   Complete  Spherical  Shell  Subjected  to  Uniformly 
Distributed/  Exponentially  Decaying  Pressure 

As  a  first  example  let  us  examine  the  case  of  a  complete 
spherical  shell  of  radius  R  and  thickness  h  subjected  to  a 
uniformly  distributed,  exponential  internal  pressure,  p^. 
Comparison  will  be  made  with  the  infinitesimal  deflection, 
two-moment  limited  interaction  (Figure  4.1)  solution  of 
Sankaranarayanan  [8.1]. 

The  deformation  field  for  this  trivial  problem  is  one  of 
uniform  radial  expansion. 

The  velocity  field  is 


w((j),t)  =  Wq  (t)  -  0  (8.1) 


There  is  no  bending  (M,   =  M^   =0);  therefore,  normality 

(J)      y 

requires  that 


N.   =  N-   =  N-  (8.2) 

cp      6      0 


where 


"o  =  "o*^ 


Dissipation  equation  (6.11)  reduces  to 
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2  2 

2  /  2Rsin(J)  N   w   d(})  =  2  /  R^sincl)  (p   -  yw  )w   dc})    (8.3) 
0  0  J       U   U 


Letting 

-t/t 


P3  =  Pq^     0 


and  evaluating  equation  (8.3),  one  finds  that 

yWQ  =  Poe"^/^0  -P^  .   (8.4) 

where  P   =  —5—  =  the  uniform  static  collapse  pressure  of  a 

C      K 

complete  spherical  shell 
and  t„  is  the  decay  time  constant. 
Integrating,  one  finds  that 


^'"O  =  -Poto^"'^''*0  -  ^c*  -^  ^0*0  '«-5) 


since 


WQ(t=0)  =  0. 


Motion  v/ill  cease  at  time  t^ ,  obtained  from  the  positive 
root  of 


-^oV'*^''''°  -  ^c^  +  Pq'^O  =  °  <S-^' 
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Integrating  equation  (8.5)  with  the  initial  condition 
WQ(t=0)  =  0  and  evaluating  at  t=t-.  yields  the  final  dis- 
placement 


t         P  t^ 

"Of  =  ^1^0*0 <^"^f''^°  +  tf  -  1'  -  ^1  -         <8-'^' 


or  through  use  of  equation  (8.6), 


"of  =  ^<^-^c'Vo  -  ^1  '^•''" 


These  results  are  in  complete  agreement  with  those  of 
reference  [8.1J. 

Section  2:   Clamped  Circular  Cylinder  Subjected  to  Uniform 
Impulsive  Load 

Consider  one  bay  of  a  rigid-plastic  cylinder  which  is 
reinforced  periodically  with  rigid  circumferential  stiffeners. 
Because  of  symmetry  the  boundary  conditions  at  each  stiffener 
are  those  of  the  clamped  condition. 

Hodge  investigated  the  response  of  such  a  shell  section 
for  infinitesimal  deflections  when  subjected  to  axisymmetric 
dynamic  pressures  [8-2]. 

Jones  considered  the  response  of  the  same  cylinder  when 
subjected  to  an  axisymmetric  impulsive  load  with  and  without 
the  effect  of  finite  transverse  deflections  [8-3] .   This 
solution  included  travelling  hinges  and  for  the  infinitesimal 
case  is  much  simpler  than  the  corresponding  dynamic  pressure 
pulse  analysis  of  Hodge  [8-2]  . 

104 


In  reducing  the  strain  rate  equations  (6.1),  equilibrium 
equations  (6.2),  and  the  energy  dissipation  rate  equation  (6.11) 
to  the  case  of  a  circular  cylinder  with  transverse  velocity 
field,  it  should  be  noted  that  ^    is  7t/2  radians,  R   is 
infinite,  and  d(J)   is  zero.   The  product  R,d<|)  is,  however, 
the  differential  longitudinal  element  dx.   Using  the 
notation 


()*=()'.  d() 


R,      dx 


and  assuming  u  to  be  zero,  the  strain  rate  equations  (6.1), 
take  the  form 

e   =  w*w*  (8.8a) 

X 

C3=|  (8.8b) 

K   =  -w**  (8.8c) 

X 

K.  =  0 


The  consistent  set  of  equilibrium  equations  (6.2)  reduces 
to 

Q   =  -M*  (8.9a) 

X       X 


(p   -  yii)  =  N*  =  0  (8.9b) 

(p.,  -  yw)  =  ^  -  (N  w*)*  -  M**  (8.9c) 

J  l\  X  X 

when  M      is  assumed  to  be  zero, 
b 

]D5 


Equation  (8.9b)  is  a  consequence  of  the  assumption 
that  u  =  0  when  d   =0. 

The  approximate  energy  dissipation  equation  (6.11) 
reduces  to  the  form 


D        N  w 

75—5-  =  E/{-| N*ww*  -  (N  w  +  M  )w**}dx 

^TTK.K       X  X       X 


-  Z  (N  w  +  M  ) .  [w*] .  =  /(p   -  yw)wdx  (8.10) 

.X      X   J  2  J 

where  the  summations  are  taken  over  i  continuous  regions  and 
over  j  hinges. 

The  velocity  profile  is  assumed  to  be  symmetric  about 
the  midspan  of  the  cylinder  (Figure  8.1). 


w(x,t)  =  WQ(t)  y;  for  0  -  X  -  L  (8.11) 


where  w^  is  the  velocity  at  x  =  L. 

Using  the  velocity  profile  (8.11),  symmetry  about  x  =  L, 
and  the  equilibrium  result  that  N*  =  0 ,  equation  (8.10)  reduces 


to 


L  N  w^x  w-  ~2^n 

2  ^  ^R^  ^^    -  2<Vx=0  -i   -     f'^x"  ^  «x'x=L  '^>    = 


L        yw  X   w  X 

2  /  (p.  -  —T-)    -f-  dx  (8.12) 

Q    3      L     L 
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From  the  two-direction  limited  interaction  yield  surface 

(Figure  4.2)  it  is  seen  that,  since  M.   =  0,  N-   =  N.  for 

V  0      0 

w  >  0. 

If  one  assumes  that  the  central  hinge  is  twice  the 
length  of  the  end  hinges  and  that  the  longitudinal  strain  is 
concentrated  equally  in  each  of  the  three  hinges,  it  can  be 
shown  that  the  ratio  of  the  average  longitudinal  strain  rate 
to  the  average  curvature  rate  in  the  hinges  is 


e   avg.      w 

^ =  -  2^  (8.13) 

K   avg. 

where  (+)  is  used  at  the  clamped  ends  and  (-)  at  x=L  when 
w  >  0. 

The  equation  of  the  x-direction  portion  of  the  yield 
surface  is 


N 
M^  =  ±Mq  [1  -  (j^)2]  (8.14) 

o 


By  normality 


dM  e 
X  _    _x 

dN  ~  1^ 

X  X 


(8.15a) 


Therefore, 


w^      2N  M^      2N  N^h 

+  0  =   +_4-0  =   ±   ^  Q^  (8.15b) 

2       ^l  ^jI   4 
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from  which 


w 

N   =  N-  ^  (8.15c) 

X     0  h 


However,  N   cannot  exceed  N-  and  the  final  result  is 
X  0 


N^  =Nq  ^  for^  ^  1  (8.16a) 

N^  =  Nq  for  j^   -  1  (8.16b) 


M^  =  ±Mq[1  -  (j^)2]        for  iT"   "  ■'■  (8.16c) 


w 

M  =  0  for  —  -  1  (8.16d) 

X  h 


where  (+)  is  used  at  the  central  hinge  when  w  >  0. 

These  results  are  the  same  as  those  of  [8-3]  for  a 
clamped  beam. 

If  equations  (8.16)  are  substituted  into  equation  (8.12) 
along  with  the  fact  that  N_   =  N.  and  the  spatial  integrations 

u  (J 

are  then  performed,  the  results,  after  some  simplification, 
are 

2M_    2M_      w^    p,    yw-  w 

and  then  become 

2M  4M      w  p  yw  w 

0+00=3 0  ^^^^      0      gj^ceeds    1.  (8.17b) 

Rh  L'^       h  2  3  h 
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#* 


It  should  be  noted  that  there  are  no  longitudinal 
membrane  stress  and  no  longitudinal  membrane  strain  rate  if 
the  transverse  deflections  are  considered  to  be  infinitesimal. 
The  kinematically  admissible  velocity  profile  (8.11)  would 
require  M   to  be  +M   and  -M^  at  the  central  hinge  and  at  the 
boundary  hinges,  respectively. 

From  equation  (8.9c)  it  can  be  verified  that  the 
generalized  stress  field 

and 


is  statically  admissible. 

The  upper  and  lower  bound  static  collapse  pressures,  P  , 
coincide. 


N     4M 
n     n 

P^  =  g^  +  ^  (8.18a) 


or 


or 


2M 
P   =  ^  (2  +  c^)  (8.18c) 


where 


2     2L^ 
c 


"  Rh  • 

P   is  the  same  as  equation  (11)  of  [8-2] 
c 
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The  dissipation  equation  (8.12)  yields  the  following 
results  for  dynamic  infinitesimal  behavior 


2Mq    2M     p     yw 


Equation  (8.19)  is  an  ordinary  linear  second  order 

differential  equation  in  w  (t) .   Equation  (8.17)  is  a  pair 

of  ordinary  second  order  differential  equations,  one  of  which 

w 
is  nonlinear  and  governs  the  motion  until  r—  reaches  unity, 

and  the  other  of  which  is  linear  and  governs  the  motion  if 

and  when  r—  exceeds  unity. 

For  the  case  of  impulsive  loading,  p^  is  considered  to 

be  zero.   The  impulse  delivered  to  the  shell  is  accounted  for 

by  the  initial  conditions  of  the  mode  approximation. 


and 


WQ(t=0)  =0  (8.20a) 


WQ(t=0)  =  Wq  (8.20b) 


where  v  is  the  Martin  and  Symonds  [8-4]  mode  factor  (defined 
below)  and  V_  is  the  characteristic  initial  velocity  deter- 
mined from  the  applied  impulse  by  the  conservation  of 
momentum. 


[  /pdt  =  yV^  ] 
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Martin  and  Symonds  showed  for  infinitesimal  deflections 
of  an  impulsively  loaded  rigid-perfectly  plastic  structure 
that  the  best  agreement  between  the  physical  deformation 
process  and  a  mode  approximation  (time-independent  velocity 
profile)  results  when  the  characteristic  mode  initial  velocity 
w   is  equal  to  a  certain  multiple  of  the  exact  characteristic 
impulsive  initial  velocity  [8-4].   It  is  assumed  herein  that 
the  same  factor  is  valid  for  finite  deflection  problems.   The 
factor,  V,  is  shown  to  be  [8-4] 


/  yu.  V.  dS 

^  =   § (8.21) 

/  yu.  u.  dS 
S    ^   ^ 


where  y  is  m.ass  per  unit  area,  u.  are  the  mode  approximation 

initial  velocity  profiles,  u.  =  V.f. (x) ,  and  V.  are  the  actual 

J  I-        '      ^  0  1  1 

impulsive  initial  velocities  V.  =  V„  g. (x) . 

For  this  particular  problem,  only  u^  and  V^  are  different 
from  zero. 


^3  =  ^0  L 


and 


^3  =  ^0 


Substituting  these  values  into  equation  (8.21),  the  mode 
factor  is  found  to  be  v  =  1.5. 
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The  motion  equations  (8.19)  and  (8.17)  can  be  rewritten 
in  the  follc-ring  form  for  the  impulsive  case 


3? 
••        c 

w-  =  -  -^^ —  for  infinitesimal  deflections 


0      2]j 


or,  for  finite  deflections 


(8.22a) 


6M   w^    3P 
^0  ^  iHJ^h^^  2^   =0  (8-22^) 

while 


and 

0    yL'    h      2y     yL^      yRh 

when 

7-—   exceeds  1  (8.22c) 

h 

where  P   is  given  by  equation  (8.18)  and  the  initial  conditions 
for  (8.22  a  and  b)  are  given  by  equations  (8.20). 

The  infinitesimal  problem  is  easily  solved  by  integrating 
equation  (8.22a)  twice  and  recognizing  that  motion  ceases 
when  w   =  0 . 

The  ratio  of  the  final  central  deflection  w...  to  the 
shell  thickness  h  is 
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h    "  leM^h  (Rh+LM  {8.23a) 


or 


^  -   3X 

h      2(2+c^)  {8.23b) 


where  X  is  the  loading  parameter 


4M  h 


Solution  of  the  finite  deflection  problem  is  more 
involved.   Solution  of  the  nonlinear  equation  (8.22b)  with 
its  initial  conditions  (8.20)  is  sufficient  so  long  as 

w 

Of  <   T 
h     "  ^* 

If  the  solution  yields 


^>  1 


then  equation  (8.22c)  must  be  solved  for  the  second  stage  of 
motion.   The  initial  conditions  for  the  second  stage  of  motion 
are  that  the  displacement  equals  h  and  that  the  velocity 
equals  the  velocity  at  the  end  of  the  first  stage  of  motion. 
Let  us  write  the  first  stage  finite  deflection  equation 
in  the  form 


Wq  +  3^  w^  +  3^  =  0  (8.24) 

6M  3P^     3M  (2+c^) 

where  B,  =  ^-j.^  2         and  3-,  =  --5 


1    viL'^h^         2    2y         yL' 
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Note  by  use  of  the  Chain  Rule  that 


0 


_   —   vv  _   -q • 

0     0  dw« 


Substitute  this  expression  into  equation  (8.24)  and 
integrate  between  the  displacement  limits  0  and  w   to  find 

^  ^^0%^  "  ^0^0^^   =  -(-^  +  32  Wq)  (8.25) 


Recalling  that 


Wq(0)  =  vVq 


and  that  w_  is  maximum  when  w^  =  0,  we  see  that 

^L  +  fi-^  ^nn,  -  Oft    =  0  (8.26) 

Om    3^    Om    2  p, 


Using  the  definitions  of  3,,   3^/  P  ,    X,  and  the  value  of  v 
for  this  particular  problem  along  with  the  formula  for  solution 
of  a  cubic  equation  [8-5] ,    it  can  be  shown  that 


?^=iF^/<¥^)'-'^'i^ 


^    [|l-/[2|Iir'    +    (|A)'    1^  (8.27) 

(First    stage    of   motion) 
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Now, 


h     h 


if  less  than  unity. 

If,  on  the  other  hand,  the  value  of  r calculated  from 

h 

equation  (8.27)  is  greater  than  unity,  find  from  equation 
(8.25)  that  the  velocity  when  r—  =  1  is 

28  h^         1 
^1  -  "^0^^^  "  Wq(t=0)  =  [v^v^  -  -^ 2&Jn]^  (8.28) 


where  t  is  time  of  second  stage  motion  (i.e.,  t=0  when  w^=h) . 

For  this  particular  problem  equation  (8.28)  can  be 
expressed  as 


9    2+3  (2+c^) 
Wp(h)  =  Wp(T  =  0)  =  Vq[|-  ^  \l^  ]  (8.29) 


The  finite  deflection,  second  stage  equation  (8.22c) 
can  be  written  in  the  form 


Wq    +   B3   Wq    +   3^    =    0  (8.30) 


where 


12Mq 

^3  =  iJiTTT 


and 


116 


6M. 


6.  = 


4    yRh 


with  initial  conditions 


Wq(t  =  0)  =  h 


and 


w  (t=0)  is  given  by  (8.29). 


The  solution  is 


w^(t)  =  :p-  sin  fx   +  (h+^)  cos  fj      -   ~ 
Of  ^3  33 


(8.31a) 


and  the  velocity  is 


w_  (t)  =  ^  -,    cos  fx 


^4 
-  f(h+^)  sin  fx 

^3 


(8.31b) 


where 


In  terms  of  this  particular  problem,  the  second  stage 
solution  can  be  expressed  as 


Wq(x) 


'3A 


(8+3c^)    .   ^    ^  4+c^      . 
-^ — ;: sm  fx   +  —. COS  fx 


(8.32a) 
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and 


^0^"^^    .  /3A    (8+3c^)       .      (4+c^)f   .   .     ,o  -,^,  , 
— j- —  =  f  /  2~  -  -^ — -p cos  fx   -  -^ 7 — - —  sm  fx    (8.32b) 


If  the  shell  does  enter  the  second  stage  of  motion,  the 
final  displacement-thickness  ratio  is 


^Of  ^   /3>^  _  (8+3cM    „.„  ^,   ^  4+c2  ^^„  ^^   _  c^ 


sin  fx,  +  -^^-^^  cos  fx,  -  i-   (8.33) 


h     /  4        6  1     4         1    4 

where  x,  is  the  first  positive  root  obtained  from  setting  the 
second  stage  velocity  equal  to  zero. 
That  is, 

^  /IX   7  (8+3c^) 


fx,  =  arctan  {  -^^^ }  (8.34) 

^  4+c2 


For  the  simply  supported  case  with  finite  deflections  it 
may  be  shown  that 


w  w     , 

N   =  2N^  ^  for  r-^  -  ^  (8.35a) 

X      0  h  h     2 

and 

w 

N   =  N^  for  c^  -  i  (8.35b) 

X     0  h     2 
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The  factor  3-.  in  equation  (8.24)  becomes 
12M. 


3,'  =  ^^  (8.36a) 


and  3~  becomes 


3M   (l+c^) 
31  =  — (8.36b) 


The  final  disolacement-thickness  ratio  becomes 


?  =  tH  v^<^>'  .  (f|)^  1^ 


+  [fl  -  /(ii^)^  +  (||)2  ]^  (8.37a) 


f    "^0  <  1 
for  j^  -  2 


or 


^0    /3X    (4  +  3c^)     .   .    ^  2+c^     ^  ^^     c^    ,„  -,_,,. 


^    "'O  >  1 
for  j^  -  2 


^  /3X    4T3c^ 


T.     ^         ^    r  V  4      12    t 
where  fx,  =  arctan  i ; 

1  o  1   2 

2+c 
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JO     X 


Equation    (8.23) 

Jones  (infinitesimal) [8-3] 

Hodge  (rectangular  pulse  lOX  static 
collapse  pressure) [8-3] 

Equations  (8.27  &  8.33) 

Jones  (finite  deflection) [8-3] 

Fig.  8.2 
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Figure  8.2  shows  the  infinitesimal  results  of 
equation  (8.23)  and  the  finite  deflection  results  of  equations 
(8.27)  and  (8.33)  for  a  clamped  cylindrical  shell  for  which 
c  =  1.   Also  plotted  are  Hodge's  infinitesimal  deflection 
results  for  a  rectangular  pressure  pulse  of  the  same  impulse 
with  peak  pressure  equal  to  ten  times  the  collapse  pressure 
and  Jone's  infinitesimal  and  finite  deflection  results  for 
impulsive  loading.   The  approximate  results  of  this  section 
are  slightly  higher  than  the  corresponding  results  from  [8-3] . 
All  results  are  for  yield  surfaces  which  circumscribe  the 
exact  yield  surface. 

It  should  be  noted  that  the  analyses  of  Hodge  and  Jones 
are  valid  only  for  c^  <  6.   There  is  no  limitation  on  c  for 
the  approximate  analysis  performed  in  this  section. 

Section  3:   Cylindrical  Panel  Subjected  to  Non-Uniformly 
Distributed  Impulsive  Load 

Reference  [8-6]  reported  some  experimental  results  on 
the  non-axisymmetric  dynamic  loading  of  cylindrical  panels. 
Of  interest  here  are  some  results  for  aluminum  6061T6  panels 
having  clamped  longitudinal  boundaries  and  free  circumferential 
boundaries  which  were  plastically  deformed  by  detonating  sheet 
explosive  which  covered  a  portion  of  the  inner  surface.   The 
midsurface  configuration  is  shown  in  Figure  8.3. 

Nominal  test  specimen  dimensions  are  shown  in  Table  8.1. 
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A 


Outline  of 
explosive  area 


Fig.  8.3 
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Table  8.1 
Nominal  Cylindrical  Panel  Dimensions  [8-6] 

L  3.0  inches 

Oq  45.0  degrees 

a  1.5  inches 

e 

8  33    degrees* 

h  0.125  inches 

0.091  inches 
0.082  inches 

'Varied  slightly  with  thickness 


The  dimensions  x.  and  6   are  mode  parameters  and  define 

0       a  ^ 

the  regions  I,  II,  and  III. 

An  approximate,  infinitesimal  deflection  solution  has 
been  obtained  using  the  general  procedure  of  Chapter  7. 

It  was  assumed  that  each  region  (I,  II,  III)  experienced 
rigid  body  motion.   Region  I  was  assumed  to  rise  vertically 
with  a  characteristic  velocity  of  w   at  x=L,  6=0. 

Region  II  was  assumed  to  rotate  rigidly  about  the  support 
line  BC  (or  AD) . 

Region  III  was  assumed  to  rotate  rigidly  about  the  line 
AB  (or  DC) . 

Because  of  symmetry  about  the  curves  x=L  and  9=0,  only 
one  quadrant  (0-x-L;  0-6-9-^)  need  be  analyzed. 
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with  the  requirement  that  w  and  w  be  continuous  at  the 
hinge  curves  between  regions,  the  displacement  rates  for  each 
region  are  found  to  be 

Region  I 

u  =  0  (8.38a) 

V  =  -w   sine  (8.38b) 
w  =  w.  cos6  (8.38c) 

Region  II 

u  =  0  (8.38d) 

V  =  n  w  {1  -  cos  (6^-6)}  (8.38e) 

a  u  u 

w  =  n  w^  sin  (6^-0)  (8.38f) 

a  0        0 


,        cose 

where  n   =       a 
a 


sinOQ-e^) 


Region  III 

Rx  w 


u  = (cose  -  cose^)  (8.38g) 

Xq  u 


x_ 
"0 


V  =  -  ^^  w^  sine  (8.38h) 

x«   0 


w  =  —  w^  cos9  (8.38i) 


The  positive  directions  of  u,  v,  and  w  are  indicated  in 
Figure  8.3. 
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Continuity  of  w  and  w  also  determines  the  equation  of 
the  hinge  curve  C^^  between  regions  II  and  III. 

sin(e  -e) 

x^^  =  x_  n   5 (8.39) 

23     0   a    cos6 

Equations  (8.38  and  8.39)  agree  with  the  results  of 
Janas  [8-7]  for  the  static  collapse  mode  velocities  and  hinge 
curve  of  a  shallow  cylindrical  panel  before  shallow  shell 
approximations  were  made. 

From  equations  (8.38)  the  in-surface  velocity  dis- 
continuities and  transverse  velocity  derivative  discontinuities 
are  found  to  be: 

Between  regions  I  and  II: 

[u]  =  0  (8.40a) 

[v]  =  (n  -  n  cos(e^-e  )  +  sine  }  w^       (8.40b) 
a    a      (J   a        a   u 

[w']  =  0  (8.40c) 


where 


and 


[w*]  =  {sine   -  n  cos(e^-e  )}  w  (8.40d) 

a    a     u   a    u 


W'   =   TT— 

3x 


w 


*  = 


9w 

ae 
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Between  regions  II  and  III: 

[u]  =  ~    (cose  -  cose.)  w_  (8.40e) 

Xq  u    u 

[v]  =   n^d  -  cos(eQ-e)  +  tane   sin(eQ-e)}  w  (8.40f) 

w 
[w']  =  -  n^  cose  (8.40g) 

[w*]  =  -  n^{cos(eQ-e)  -  tane   sinCe^-e)}  w    (8.40h) 


Between  regions  I  and  III: 

*     R  • 

[u]  =  — -  (cose  -  cose^)  w^  (8.401) 

Xq  u    u 

[v]  =  0  (8.40J) 

[w']  = cose  (8.40k) 

^0 

[w*]  =  0  (8.401) 


At  the  clamped  boundary: 

[u]  =  [v]  =  [w']  =  0  (8.40m) 


[w*]  =  -  ti^Wq  (8.40n) 


If  the  curvature  rates  for  a  cylinder  are  found  from 
equations  (6.1)  and  then  integrated  across  a  hinge,  the 
relative  rotation  rate  components  can  be  shown  to  be: 
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and 


fi^  =  [W]  (8.41a) 


.^  =  tv]-[w*]  (g^^^^j 


which  may  be  evaluated  using  equations  (8.40). 

With  the  six-dimensional  approximate  yield  surface 
described  in  Chapter  7  (Figure  7.2)  normality  considerations 
for  w^>0  require  the  following  generalized  stress  values. 

Between  regions  I  and  II: 


M   =N   =M^=N^=0  (8.42a) 

X   X   xe    xe 


Between  regions  I  and  III 


M   =  M_;  N   =  N^ 
X     Ox     0 


"e  =  Ng  =  M^g  =  "xe  =  °  <^-^2b) 


Between  regions  II  and  III 


M   =  M„  =  2M  ^  =  M_ 

X    e    x0    0 


N   =  N^  =  2N  .  =  N^  (8.42c) 

X     9      x6     0 


At  the  clamped  boundary: 

M   =N   =M„=N„=N^=0  (8.42d) 

X     X     x9     x9     9 
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If  the  appropriate  deformation  rates  (8.40  and  8.41) 
and  generalized  stress  values  (8.42)  are  substituted  into 
the  internal  energy  dissipation  rate  expression  (equation 
7.6)  and  the  integration  is  carried  out,  it  may  be  shown  that 
the  dissipation  rate  for  one  quadrant  of  the  panel  divided 
by  the  characteristic  velocity,  w  ,  and  the  yield  moment, 
Mq  is 


D.  .  ri 

-^^  =   (l-Xq)   {^(n^[i  -  cos^Q-e^)]   +  sine^)   +  —■} 


-0^0 


,2  T,  LTI 


and 


"■  §T  ^^i^^o  -  ^0  ^°^^o^   ^  f^  ^^^^0  ^  -^ 


+  -r--^  ne.-e^)  -  cos^e^  (tane_-tan0^)}  +  ^  ne.-e^) 

h  Oa  0  0.  azOa 


sece„+tan9^  x^ 


4x 

+   -vr-     n^    cose^    { (tan0_-tane    )    -    ?    cosi 
n  a  u  u  a 


sec6_+tan9_ 

[Zni ^^4-r s^)    +   tane,.      sece_    -    tanB^      sece    ]  }  (8.43) 

sece  +tan6         0       0       a      a 


Let  us  call  the  right  hand  side  of  this  expression  F. 


F  =  F(L,  Xq,  Gq,  9^)  (8.44a) 


D.  =  F  M^  w^  (8. 44b) 

1       0   0 
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The  external  energy  dissipation  rate  (equation  5.1) 
becomes 

D      =    /    [  (p--yw)w   -    yiiu   -    yvv]    dS  (8.45) 

S   quadrant 

The  indicated  integration  must  be  performed  in  each 

region,  and  the  expressions  for  the  pressure  and  velocities 

within  each  region  depend  on  the  assumed  size  of  the  central 

region  (Region  I)  relative  to  the  loaded  area  (defined  by 

I      and  e  ) . 
e      e 

Since  u  and  v  are  defined  in  terms  of  w- ,  equation  (8.45) 
can  be  expressed  as 

Dg  =  (P3g  -  vewQ)wQ  (8.46) 


where 


g  =  g(L,  Xq,  l^,    Bq,  6^,  0^) 


and 


e  =  e(L,  Xq,  £^,  6^,  &^,  5^) 


The  governing  equation  of  motion  is  found  by  equating 
the  internal  and  external  energy  dissipation  rates  and  has 
the  form 
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FMq  =  P3g  -  yeWQ  (8.47) 


The  upper  bound  collapse  pressure  is 


c    g 


v/hich  for  a  given  panel  with  fixed  loading  geometry  varies 

with  the  mode  shape  parameters  x-  and  6  . 

^  0      a 

The  pressure  p..  is  zero  everywhere  outside  of  the  loading 
area  (Figure  8.3)  and  g  is  found  by  evaluating  the  integral 


g=/wdS  (8.48) 

S        ^ 
e 

where  S   indicates  the  loaded  portion  of  the  midsurface  area, 
e  ^ 

Equation  (8.48)  consists  of  contributions  from  each 

region,  and  the  limits  of  integration  within  each  region 

depend  on  the  parameters  x_  and  9  .   The  form  of  the  equation 

for  the  values  of  x_  and  6   indicated  in  Figure  8.3,  for 

0       a  ^  ' 

example,  would  be 

^a      L  ,  ^e      L 

g=/         /         Rwdxde-i-/         /  R  w.^.^    dx   de 

0         ^0  ^a      ^23 


^e  ^23 

+      {/      /        +      /  /      }         R  Wjjj   dx   de  (8.49) 

6  x- 
a        0 


e 

a 

L 

{/ 

/ 

0 

X 
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where  the  subscript  I  indicates  region  I,  etcetera.   Equation 

(8.39)  is  used  to  define  x_^  in  the  limits  of  integration. 

The  optimum  values  of  the  mode  parameters  are  those 

which  give  the  lowest  static  collapse  pressure.   Because  of 

the  complicated  forms  of  F  and  g,  parameter  optimization  was 

accomplished  by  computing  the  collapse  pressure  for  discreet 

combinations  of  x_  and  0  .   In  all  cases  it  was  found  that 

0       a 

the  minimum  collapse  pressure  occurred  for  x-=L  and  for  6 

^    ^  0  a 

somewhat  less  than  6  .   Region  I  degenerated  into  a  curve 
as  shown  in  Figure  8.4. 

Now  e  may  be  determined  as 

-j        a  L       0    23     ^     ^     ^ 

e  =  .  •>    [  /   /   +   /    /    ]  (uu  +  vv  +  w\v)^_^  R  dx  de 
"'O^'O     0   0      0    0  ^^^ 


n    ®0  L    .    .    . 

+  .  ..    /  /    (uii  +  vv  +  ww)^T-  R  dx  d9  (8.50) 

^0^0   0  X,,  ^^ 

a    23 


where  the  mode  velocities  and  accelerations  for  the  first 
set  of  integrals  are  those  of  region  III  and  those  of  the 
final  integral  are  of  region  II.   Evaluation  of  this 
expression  leads  to 


6    2  e 

1^  =  ~  +  ^(^  +  k   sin20   -  2cose^  sin0   +  0   003^0^) 
RL3L'^2^       a         0      aa       0 
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Fig.     8.4 
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n    R^  sine 

+   -^rr-    [7  (e_-e^)sine_    +   — j-^    (sin2e-    -    sin2e    ) 
■lj  u      a  u  4  u  a. 


•sin2eQ    (sineQ    -    sinS^)    +    Og'^a^    sinBQ    cos^eQ 


cos6 


+ — ^ —  (cos2e   -  cos2e  )  -  2cos^e   (cosBq  -  cose  ) 


Si 


+  sinS^  cos^e^  (tane.  -  tane   -  6-  +  e  ) 
0      0      0       a    0    a 


seco-. 

-2(6.-0  )  sine.  +  2cose_  Zn 5-^+  2sin2e.  (sine^- 

0   a      0        0     sece  0       0 


sine  )  +  2cos26_  (cose«  -  cose  )  -  sin2e-  • 
a  0       0       a         0 

sece„+tane     cos^e  sece- 
rn  E^ a^ o [^(tan^e-  -  tan^e   -  £n A} 

sece  +tane       3     lz  v      q        ^  sece 


(8.51) 


From  the  definition  of  the  Martin  and  Symonds  mode  factor, 
V,  equation  (8.21)  it  can  be  seen  that  for  constant  density, 
\i ,    the  denominator  of  v  is  yeV^.   The  impulsive  velocity,  V^, 
is  equal  to  the  explosive  impulse  divided  by  the  mass  of  the 
panel  portion  included  within  a  radial  projection  of  the  sheet 
explosive  boundaries. 
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For  the  degenerate  mode  shape  shown  in  Figure  8.4/  the 
numerator  of  v  is 

e         e   x_- 

a  e    23 

yJV^  =  [  /   /     +     f        f  ] (yv^     V. )  Rdx  de 

0    L-£       e    L-£        III 
e      a     e 

+  /    /     tjV      V.  Rdx  de  (8.52) 


where  the  V.    are  the  mode  velocity  components  in  region  A 
and  the  V.  are  the  actual  impulsive  velocity  components. 
Performing  the  indicated  integration,  one  finds  that 


J  =  1^  {L^e^  -  (L-£^)2  6^}  +  RLn^  { (e^-e^)  sine, 

2L       a       e     e       a    e   a      i 


+  cose  -  cose  -  T]   [^  (sin^e^-cos^e,,)  (e  -e  ) 

e       a     a  '^       0      0    e   a 


+  sine^  (cose   -  cose  )  +  cose^  (sine   -  sine  ) 
0      e       a        0      a       e 

sece   +  tane  cos^e,. 

+  cose^  £n Q^  ^  ^ ~  +  ~    (tane  -tane^)]}     (8.53) 

0    sece   +  tane        2         e      a 
a        a 


For  the  impulsive  solution  p^  is  assumed  zero  and  the 
initial  characteristic  velocity  is  assumed  to  be 


WQ(t=0)  =  vVq 


,  J 

where   v  =  — 
e 
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The  equation  of  motion  (8.54)  becomes 


™0 
"O  =  -W-  (8.55) 

"O  =  -i:¥-  *  *   ^^0  '8-56) 


Motion  ceases  at  t^  when 


"o'*f'  =  ° 


and 

yevV 


tf  =  -Fm7  '8-") 


The  final  displacement-thickness  ratio  is  found  by 
evaluating  (8.58)  at  t^. 


^Of    ^^'   ^0 


2eFM  h 


(8.59) 


Figure  8.5  shows  the  final  central  permanent  displacement- 
thickness  ratios  observed  during  the  experiments  on  the  four- 
teen aluminum  panels  reported  in  [8-6]  plotted  as  a  function 
of  a  non-dimensional  impulse  parameter,  i|i . 

yR  vg 

ij;  E   -   ^  (8.60) 

M 
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The  same  figure  contains  a  plot  of  points  calculated  from 
equation  (8.59)  for  panels  of  the  three  different  nominal 
thicknesses  reported  in  [8.6].   Within  plotting  accuracy,  all 
calculated  points  fell  on  the  same  straight  line.   The  lower 
line  of  calculated  results  is  based  on  the  approximate  yield 
surface  which  circumscribes  the  exact  six-dimensional  Tresca 
yield  surface.   The  upper  line  of  calculated  results  is  based 
on  the  approximate  yield  surface  whose  scale  is  0.618  that  of 
the  circumscribing  approximate  yield  surface.   The  smaller 
yield  surface  may  not  inscribe  the  exact  yield  surface,  but, 
as  was  noted  in  Chapter  7,  no  approximate  yield  surface  which 
reduces  to  the  two-moment  limited  interaction  yield  surface 
and  has  a  scale  factor  greater  than  0.618  can  possibly  inscribe 
the  Tresca  yield  surface. 

The  fact  that  the  smaller  yield  surface  gave  better 
results  for  this  problem  than  the  larger  yield  surface  is 
not  surprising.   A  large  portion  of  the  hinges  has  generalized 
stresses  corresponding  to  the  corners  marked  A  and  B  in 
Figure  7.2.   Corners  of  a  circumscribed  yield  surface  lie 
outside  of  the  exact  yield  surface,  while  the  corners  of  the 
inscribed  yield  surface  may  lie  on  the  exact  yield  surface. 

Section  4:   Spherical  Cap  Subjected  to  Uniformly  Distributed 
Impulsive  Load 
As  a  final  example  let  us  consider  the  case  of  a  deep 
spherical  cap  with  clamped  boundary  subjected  to  a  transverse 
dynamic  pressure  load. 
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Giannotti  [8-8]  reported  the  results  of  some  experiments 
on  clamped  spherical  caps  which  were  loaded  by  detonating 
sheet  explosive  which  covered  most  of  the  internal  surface. 
Of  interest  here  are  the  results  of  some  experiments  on 
6061T6  aluminum  caps.   Nominal  dimensions  are  listed  in 
Table  8.2  and  indicated  in  Figure  8.6. 

Table  8.2 
Nominal  Spherical  Cap  Dimensions  [8-8] 


*0 

*e 

R 

h 

(degrees) 

(d( 

agrees) 

(inches) 

(inches) 

90 

89 

2.58 

0.115 

90 

89 

2.58 

0.176 

90 

89 

2.58 

0.237 

60 


59 


2.58 


0.117 


45 


44 


2.58 


0.119 


Onat  and  Prager  considered  the  static  collapse  problem 
of  a  deep  clamped  shell  subjected  to  a  uniformly  distributed 
pressure  [8-9].   With  an  assumed  velocity  profile  of  the  form 


WQ((|),t) 


u  =  0 


W   (t)   (cose})  -  COS(})q) 

(1  -  coscf)^)     ^ 
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(8.61a) 
(8.61b) 


Fig.     8.6 
139 


where  w^  is  the  transverse  velocity  at  i>=0 ,    and  as  a  conse- 
quence of  using  the  Tresca  yield  surface,  the  state  of 
generalized  stress  was  found  to  be  membrane  (i.e.,  M.  = 
M„  =  0)  except  in  a  narrow  region  near  the  boundary, 
(J)*  -  c|)  -  (})  .   The  angle  (p*    is  defined  by  the  relation 

2cos(J)^ 
cos(})*  =  r-^  ,  (8.62) 

2-1 


derived  from  the  norinality  requirements  of  plasticity. 

The  upper  and  lower  bound  collapse  pressures  obtained 
by  Onat  and  Prager  for  shells  in  the  range  of  cf)^  and  r- 
listed  in  Table  8.2  agreed  within  15  percent  or  better. 

The  loading  used  by  Giannotti  [8-8]  appeared  sufficiently 
similar  to  uniform  loading  to  use  the  Onat  and  Prager  collapse 
profile  (equation  8.61). 

For  comparison  purposes  the  spherical  cap  problem  was 
solved  for  infinitesimal  deflections  and  for  finite  deflections 
Two  infinitesimal  solutions  with  the  Onat  and  Prager  velocity 
profile  were  obtained,  one  using  the  Tresca  yield  surface  and 
one  using  the  two-moment  limited  interaction  approximate  yield 
surface  (Figure  4.1).   Lower  bounds  on  the  final  displacement 
were  also  obtained  using  the  method  of  Morales  and  Nevill 

[8-11]  for  a  perfectly  plastic  material. 

First  let  us  examine  the  finite  deflection  solution. 

In  order  to  simplify  the  finite  deflection  analysis,  the 

two-direction  limited  interaction  yield  surface  (Figure  4.2) 
was  utilized  in  place  of  the  exact  Tresca  yield  surface. 
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The  strain  rate  and  curvature  rate  equations  (6.1)  for 
a  spherical  cap  with  the  velocity  field  of  equation  (8.61) 
reduce  to 

Wq  w  sin^cf) 

^A  -   5Tt — ^^^A  \  t^cos(})  -  cos(}).  +  ■= T— }  (8.63a) 

(p    R(l-cos(})-)      ^       ^0    l-cos({)- 

W   (COSCJ)  -  COS(})-) 

^e  =   r(i-cos^q)  (8-^3b) 

Wq  COS(j) 

K,  =  Kq  =  (8.63c) 

'^  ^        r2  (l-cos(t)Q) 

The  last  term  of  equation  (8.63a)  is  the  only  finite 
deflection  contribution. 

If  we  assume  that  the  two  hinges'  rotation  rates  can  be 
represented  by  average  curvature  rates,  each  one  averaged  over 
one  half  of  the  arc  RAc})*,  in  conjunction  with  a  meridional 
strain  rate  averaged  over  the  arc  RA4)*,  it  can  be  shown  from 
the  normality  requirements  of  plasticity  that 

N,      _,   sind)^  -  sind)*  -  Ad)*  cosA. 
/  (p  .   _  R   0 0 

^N^A  ~  h  sin<})Q 


^0 


A4,*    sin2^Q  -  sin2(j)* 


+  _ii  — ± . '. (8.64a) 

sin(})Q  (l-cos(t)Q) 


and 
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N,      _,   sin(})-  -  sinft)*  -  Ad)*  cosd),, 
,_9\   _  R      0 ]_0 

'^Q   B        h  sincf)* 


Ac},*    sin2(})Q  -  sin2(})* 


+  rr-     ^  .  ,^  ,-, ^-j— X (8.64b) 

h        sine})*  (l-cos(})^) 


where  A  and  B  refer  to  hinges  A  and  B  of  Figure  8.7. 

The  bending  moment  associated  with  each  hinge  may  be 
found  from  the  yield  surface  equations 


M,  N^ 


_i  =   ±(1 ^y  (8.64c) 

0  N; 


where  (-)  is  used  for  hinge  A  and  (+)  for  hinge  B  when  w_>0. 
For  (})-(})* 


N^  =  N^  (8.64d) 

and  • 

M,  =  0  (8.64e) 


Circumferential  bending  is  important  only  in  the  same 
region  as  meridional  bending  and  it  can  be  shown  that  for 
Wq>0. 


-1  =  £^  (1 -0)  (8.65a) 

N     h        cos(}) 

and 
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^0*  =  0„-0 


^ 


Fig.    8.7 
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^e  _  ,  _  4r2  ,,   ^°^*o.2 


JJ-  =   ^   -   -^      (1  -  r^)  (8.65b) 

M-.        h        cos())  v«.v^-'x^/ 


0 


for  (})*  -  (j)  -  (()^ 


while 


Nq  =  Nq  (8.66a) 


and 

M„  =  0  (8.66b) 


for  (})-(})* 


If  the  velocity  profile  (equation  8.61),  pressure,  and 
the  generalized  stress  values  (equations  8.65  and  8.66)  are 
substituted  into  the  external  dissipation  rate  expression 
for  shells  of  revolution  (equation  6.11),  and  the  indicated 
integration  performed,  it  can  be  shown  that 

Re,  Rc^ 

Re,  ^^1^2  ^^1^2 

+  P  {a_+a^  T-TT X 2a^  .  ,7  -y^  -  2a-  — r }   w 

c   2   5  hd-y)      6  hd-y)       7    h       0 

^2  ^2  ^2 

c    5  h(l-y)     6  h(l-y)        7  h'^    0 


^8^3  "  ^9^^0  (8.67) 
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where 


Rh 


P   is  the  Onat  and  Prager  lower  bound  collapse  pressure 


_  sin(j)' 


^    sin(|)Q 


sine})-.  -  sine})*  -  Ac})*  cose}), 
c,  = 


'1  sin(})^ 


^A*    sin2(I)Q  -  sin2(})* 


2 


2        sincf)_  (l-cos(}>_) 


l-cos^d)*  ,     ,      ,+       ■ 

a,  =  = ^^—  +  cos(|)_  coscj)*  -  cos(()^ 

=    1  .1-cos^cD*   ^  cos^q(cos(D*-1) 

^2    R(1-cos(|)q)  ^32         ^ 


a-  =  r-  [i'  (3C0S^(})-  +  COS^C})*)   -  2C0S(})    COScJ)*  +  COS^(j)_ 


£n^^] 

COS(})q 


^4  "  3^[i  (cos2(|)*  -  cos2({)q)  -  IP  a^] 


sine})*  (sin(})_  -  sine})*)  (coscJ)*  -  cos({)_) 
^5  "  2R(1-cosc{)q) 


^6  ~  8R  ^^"^*  (sin(})Q  -  sine})*) 
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h  sin^cf). 
a_  = 


7       8R 


ag  =  -^  (1-cos(})q)^ 


ag  =  ^   (1-cos(})q)  3 


and 

y  is  the  mass  density  per  unit  midsurface  area  . 

By  defining  new  constants,  equation  (8.67)  can  be 
written  in  the  form 


where 


a^b^  is  the  coefficient  of  P  , 
9  1  c 


a„b^  is  the  coefficient  of  P  w. , 
3  2  C  0 


and   a^b^  is  the  coefficient  of  P  w^  on  the  left  hand  side  of 
9  3  c  0 

equation  (8.67) . 


^4  =  ^ 
4    ag 


Equation  (8.68)  has  been  solved  for  the  impulsive  loading 
case  by  the  method  of  successive  approximations  [8-10]. 


146 


The  pressure  p-  is  assumed  to  be  zero  and  the  initial 
characteristic  velocity  is 


WQ(t=0)  =  vVq  (8.69) 


For  the  partial  loading  of  the  spherical  cap  in 
conjunction  with  the  assumed  velocity  profile  (equation  8.61) 
the  Martin  and  Symonds  mode  factor  becomes 


COSct)    -  COS(J)_ 


The  actual  impulsive  initial  velocity,  V- ,  is  equal  to 
the  explosive  impulse  divided  by  the  mass  of  the  spherical 
cap  included  within  a  radial  projection  of  the  sheet  explosive 
boundary. 

Using  the  second  approximation  to  the  solution  of 
equation  (8.68)  one  finds 


w^ (t)  =  J,  sinAt  +  J^  cosAt  +  J^  sin^At  +  J-  cos^At 


+  J^  sin2At  +  J^  t  sinAt  +  J_  t  cosAt  +  Jq  (8.71) 

5  6  /  o 


where 


A  =    /PbT 
c    2 


vV  b    b 

1  A  ,2 

^2 
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and 


^    ^2    3b^      3b^A2 


b^    v^V      2b ^ 
^    ^"2    A^       b| 


b     2v^Vg    bj 


J  ^  .  t'li's^^o 


^     3b^A 


''i'^3^ 

's—^ 


b^b3VVp 


^1    '^I'^S 


The  velocity,  w  (t) ,  can  be  found  by  differentiating 
equation  (8.71)  with  respect  to  time.   Motion  ceases  at  the 
time,  t^,  when  w.  =  0 .   The  time  t^  is  thus  the  first  positive 
root  of  the  transcendental  equation  obtained  by  setting  w. 
equal  to  zero. 

The  final  central  displacement,  w  _,  is  obtained  by 
evaluating  equation  (8.71)  at  time  t,;. 
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Computed  results  show  less  than  one  percent  difference 
between  the  first  and  second  approximations. 

These  approximate  analytic  results  are  themselves  quite 
complicated,  due  primarily  to  the  presence  of  finite  deflec- 
tion terms.   In  order  to  determine  whether  or  not  the  additional 
complexity  is  justified,  the  problem  was  also  solved  for 
infinitesimal  deflections  alone. 

If  the  same  velocity  field  and  the  exact  Tresca  yield 
surface  are  assumed,  and  if  a  hinge  is  allowed  at  (j)»,  the 
dissipation  rate  expression  for  shells  of  revolution 
(equation  6.11)  produces  the  following  infinitesimal  results 


P^(a^  +  2a3  +  2a^  +  a^)  =  aQP3  -  a^w^  (8.72) 


Equation  (8.72)  is  a  linear,  ordinary  second  order 
differential  equation. 

For  the  impulsive  case,  p^  =  0.   Integrating  twice  with 
respect  to  time  and  using  the  initial  conditions  that 


w^(t=0)  =  vVq 
and 


WQ(t=0)  =  0, 
one  finds  that 


P  (a,+2a^  +  2a.+a_,)  t^ 

w  (t)  =  -   ^   ^  ^J ^— ^ +   vV_t  (8.73) 

o  2a-,  0 
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Finding  t^,  from  the  fact  that  motion  ceases  when 
w_(t^)=0  and  substituting  into  (8.73)  leads  to  the  final 
permanent  central  displacement-thickness  ratio 


^Of  _        ^9  ^'  ^g  ,,  ^,, 

h      2P  h{aT+2a^+2a^+a^)  la./^^ 

c   1    3    4   7 


If  the  Onat  and  Prager  velocity  profile  is  used  in 
conjunction  with  the  two-moment  limited  interaction  yield 
surface  the  following  infinitesimal  deflection  results  are 
found 

d-cosc}).)^    ,  p,    yw 

P^  { ^-^  +  %    {l-cos  =  *g)}  =  (/  -  -^Oxi-cos*,)^ 


(8.75) 


The  impulsive  final  central  displacement-thickness  ratio 


is 


w-^     y  v^    V^ 

Of  = 0 (8.76) 

^     3P  hn+^  l+cos(j)o 
c     2R  i-cos(})Q^ 

Morales  and  Nevill  [8-11]  have  developed  a  lower  bound 
displacement  technique  for  rigid-plastic  structures  subjected 
to  impulsive  loading.   The  method  is  based  on  the  use  of  a 
kinematically  admissible  velocity  field. 

The  lower  displacement  bound  obtained  from  the  Onat  and 
Prager  velocity  field  and  the  Tresca  yield  surface  is 
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«^ 


^   >   ^9  ^'  ^0 
h      3P  h(aT+2a-+2a,+a^)  (8.77) 

c    1    3    4   7 

The  lower  bound  is  seen  to  be  two-thirds  the  value 
predicted  by  the  Tresca  yield  surface  infinitesimal  solution. 

Comparison  of  experimental  and  approximate  theoretical' 
results  are  shown  in  Figures  8.8  through  8.13.   These  figures 
show  final  permanent  central  displacement-thickness  plotted 
as  a  function  of  impulse  parameter  X. 

yR^V^     4pR2v2 
^  =  "TyT  "   a^h^  (8.78) 

where  p  is  density  based  on  volume. 

Figure  8.13  is  a  combined  plot  of  the  information  from 
Figures  8.8,  8.9,  and  8.10.   This  figure  shows  that 
if  one  combines  the  experimental  results  for  the 
individual  series  of  90  degree  caps  the  complete  set  of 
results  approximates  a  curve. 

Each  figure  (except  8.13)  has  a  plot  of  experimental 
results,  Tresca  infinitesimal  results,  two-moment  limited 
interaction  infinitesimal  results,  Tresca  lower  displacement 
bound,  and  finite  deflection  results. 

It  is  immediately  obvious  that  in  the  range  of  observed 
displacement-thickness  ratios  the  effect  of  finite  deflections 
on  the  analytic  solution  is  negligible. 
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It  is  also  obvious  that  agreement  with  experimental 
results  is  only  fair.   Particularly  noteworthy  is  the  fact 
that  some  experimental  results  for  the  60  degree  caps  are 
below  the  infinitesimal  deflection  lower  bound  for  a  perfectly 
plastic  material  and  all  experimental  results  for  the  45  degree 
cap  are  below  the  lower  bound. 

The  90  degree  cap  approximate  theoretical  results  can  be 
plotted  in  a  manner  which  show  final  displacement  to  be 
essentially  independent  of  thickness  for  a  given  initial 
velocity.   This  behavior  is  indicative  of  membrane  behavior. 
Indeed  bending  does  have  very  little  influence  in  the  analysis 
of  this  problem.   The  experimental  results  do  not  show  this 
independence  of  thickness,  which  indicates  that  bending  may 
be  of  greater  importance  than  assumed. 

In  the  interests  of  simple  analysis  it  was  assumed  that 
u=0.  In  order  to  include  greater  bending  dependence,  it  may 
be  necessary  to  utilize  a  more  complicated  velocity  field. 

The  simple  infinitesimal  results  obtained  for  this 
problem  may,  however,  serve  as  useful  engineering  estimates 
of  shell  response. 
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CHAPTER  9 
SUMMARY  AND  CONCLUSIONS 

Exact  solutions  of  the  response  of  shells  dynamically 
loaded  above  the  static  plastic  collapse  load  are  quite 
complicated.   General  approximate  theoretical  procedures  are 
developed  herein  for  arbitrarily  shaped  shells  with  arbitrary 
boundary  conditions  and  with  arbitrary  loading  made  of  rigid- 
perfectly  plastic  material.   The  influence  of  finite  trans- 
verse deflections  is  considered. 

Approximate  results  were  compared  with  exact  analytic 
solutions  for  complete  spheres  and  circular  cylinders  and 
also  with  experimental  results  for  cylindrical  panels  and 
spherical  caps. 

The  approximate  results  agreed  exactly  with  the  exact 
analysis  for  a  complete  sphere  and  agreed  quite  well  with 
the  exact  analysis  of  a  complete  cylinder  loaded  axisymmet- 
rically.   The  approximate  results  for  partially  loaded 
cylindrical  panels  also  agreed  favorably  with  experimental 
results. 

The  approximate  results  for  explosive  loading  of 
spherical  caps  underpredicted  some  experimental  results  and 
overpredicted  others.   The  variation  v;as  as  great  as  60%  of 
the  predicted  results. 

More  scatter  was  observed  in  the  experimental  results 
for  spherical  caps  than  has  been  observed  in  related  dynamic 
experiments  (e.g.  the  cylindrical  panel  results).   Some  of 
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the  spherical  cap  experimental  results  are  below  a  lower 
bound  on  displacement  for  infinitesimal  deflection  of  a 
perfectly  plastic  material. 

It  appears  from  the  approximate  results  obtained  that 
the  influence  of  finite  deflections  is  quite  small  if  the 
shell  geometry  requires  the  development  of  sizeable  membrane 
stresses  with  infinitesimal  deflections.   Therefore,  it  seems 
reasonable  to  use  infinitesimal  theory  to  predict  permanent 
displacements  for  such  shells. 

It  would  be  desirable  to  develop  the  six-dimensional 
generalized  stress  space  Tresca  yield  surface  at  least  to 
determine  the  size  ratio  between  inscribed  and  circumscribed 
approximate  yield  surfaces. 

In  future  experiments  concerning  the  dynamic  behavior  of 
explosively  loaded  shells,  it  would  be  instructive  to  deter- 
mine photographically  whether  or  not  the  assumed  impulsive 
initial  velocity  profile  is  correct. 

Finally,  it  would  be  worthwhile  to  apply  the  techniques 
of  this  thesis  to  other  shell  geometries  and  conditions  as 
more  experimental  results  become  available. 
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APPENDIX  A 

RELATION  BETWEEN  THE  TWO-DIRECTION  LIMITED  INTERACTION 

YIELD  SURFACE  AND  THE  TRESCA  YIELD  SURFACE 

In  Chapter  4  it  was  stated  that  the  two-moment  limited 
interaction  yield  surface  and  the  two-direction  limited 
interaction  yield  surface  (Figures  4.1  and  4.2)  circumscribe 
the  exact  Tresca  yield  surface  and  that  if  each  were  reduced 
in  size  to  0.618  times  original  size  they  would  inscribe  the 
Tresca  yield  surface. 

This  appendix  contains  as  a  sample  calculation  the 
limiting  condition  for  the  two-direction  limited  interaction 
yield  surface. 

The  procedure  used  to  inscribe  an  approximate  yield 
surface  within  the  Tresca  yield  surface  is  outlined  in 
reference  [A-1] .   It  consists  of  assuming  a  state  of  stress 
corresponding  to  a  vertex  of  the  approximate  yield  surface 
and  then  finding  the  common  positive  factor  by  which  the 
stresses  must  be  multiplied  in  order  for  the  vertex  to  lie  on 
the  Tresca  yield  surface.   The  smallest  value  of  the  scale 
factor  must  be  used  to  inscribe  the  approximate  yield  surface 
within  the  Tresca  yield  surface. 

Since  yield  surfaces  exhibit  symmetry,  not  every  combina- 
tion of  vertices  has  to  be  examined.   Symmetry  of  the  two- 
direction  limited  interaction  yield  surface  is  such  (see 
Figure  4.2)  that  only  the  first  quadrant  of  each  curve  need 
be  examined. 
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A  limiting  vertex  of  the  approximate  surface  is  the 
state  of  stress 


"2  =  N^  =  1 


1  Mq 


M 
2 

2    Mq 


If  these  values  are  multiplied  by  the  common  scale 
factor  E,    and  the  values  n,  =  m  =  0  and  n„  =  m,  =   E,    are 
substituted  into  the  defining  equations  of  the  Tresca  yield 
surface  [A-1] ,  the  limiting  equation  is  found  to  be 


m   -  m   =  1  -  (n^-n  ) ^  (A.l) 


which  becomes 


C  =  1  -  C^  (A. 2) 


Solving  for  the  positive  value  of  the  scale  factor  yields 


-1  +  /5 
^  =  -^        =  0.618  (A. 3) 
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